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The methods of mode decomposition and Fourier analysis of classical and quantum 
fields on curved spacetimes previously available mainly for the scalar field on Friedman- 
Robertson- Walker (FRW) spacetimes are extended to arbitrary vector bundle fields on 
general spatially homogeneous spacetimes. This is done by developing a rigorous unified 
framework which incorporates mode decomposition, harmonic analysis and Fourier anal- 
ysis. The limits of applicability and uniqueness of mode decomposition by separation of 
the time variable in the field equation are found. It is shown how mode decomposition 
can be naturally extended to weak solutions of the field equation under some analytical 
assumptions. It is further shown that these assumptions can always be fulfilled if the 
vector bundle under consideration is analytic. The propagator of the field equation is 
explicitly mode decomposed. A short survey on the geometry of the models considered in 
mathematical cosmology is given and it is concluded that practically all of them can be 
represented by a semidirect homogeneous vector bundle. Abstract harmonic analytical 
Fourier transform is introduced in semidirect homogeneous spaces and it is explained 
how it can be related to the spectral Fourier transform. The general form of invariant 
bi-distributions on semidirect homogeneous spaces is found in the Fourier space which 
generalizes earlier results for the homogeneous states of the scalar field on FRW space- 
times. 
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1. Introduction 

As long as mankind is not in possession of a successful and commonly accepted 
quantum theory of gravitation (and possibly even thereafter), the quantum field 
theory on curved spacetimes (QFT in CST) is an adequate and consistent theo- 
retical framework for astronomy and cosmology. Vicinities of black holes and the 
early epoch of the universe are two prominent physical situations where gravity 
is sufficiently strong so that its influence on the quantum field theory cannot be 
neglected. At the same time, in these situations the gravity is sufficiently uniform 
(i.e., its local fluctuations are negligible) to be considered classical and interact- 
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ing with matter quantum fields mainly macroscopically. This semiclassical picture 
is captured by letting quantum fields propagate on a curved spacetime. The back 
reaction of the matter on gravity is described by the semiclassical Einstein equa- 
tion, where gravity feels only the expectation values of quantum entities. This is 
the domain of quantum field theory on curved spacetimes, apart from its intrinsic 
appeal for the beauty and variety of fundamental problems it poses in mathemati- 
cal physics. QFT in CST adopts the more modern algebraic quantum field theory 
setup and, as appropriate to mathematical physics, attempts to be as axiomatic 
and deductive as possible and mathematically rigorous. The disadvantage of such 
an approach is the extreme difficulty of producing explicit ready-made results which 
can be applied on the observational level, and each such result can be considered 
as a remarkable success. For instance, it was not until 2010 when the authors of [1] 
obtained the first completely rigorous and at the same time explicit description of 
cosmological particle creation in states of low energy. This was done under several 
assumptions which can or cannot be considered realistic in cosmology. Namely, the 
Klein-Gordon field on FRW spacetimes was chosen, and homogeneity and sufficient 
regularity was stipulated a priori for the desired state of low energy. When one tries 
to step a bit beyond these restrictions one immediately faces severe mathematical 
difficulties along the entire way from the very setup until the final expressions. The 
reason is that the chain of results used in these constructions has been obtained 
only under the above mentioned assumptions. The aim of the current work is the 
extension of some of those mathematical methods to a generality where they can 
be applied for practically all realistic cosmological situations. To which extent this 
program has been successful will become clear below. 

A primary tool for obtaining explicit constructions are geometric symmetries. 
After publishing his eminent work on general relativity, Einstein declared he had 
no hope to see explicit solutions of his equation in the near future. It was the 
rich symmetry of the FRW spacetimes that allowed Friedmann to find first explicit 
solutions and thus to dispel the despair of Einstein shortly after his publication. 
This instant can be considered as the birth of modern mathematical cosmology, 
which until today remains one of the main appliers of explicit solutions in general 
relativity. One of the merits of geometric symmetries is the possibility of the sep- 
aration of variables in the field equation which helps to obtain explicit solutions. 
The mode decomposition of the solutions of the field equation (also referred to as 
the Fourier method in PDE, or expansion into harmonic oscillators in physics) was 
probably first applied in the cosmological context by Parker [2] who performed it 
on the flat FRW spacetime. The idea of the method is that one tries to separate the 
time variable in the field equation, and looks for solutions as linear combinations 
of products X{x)T{t) where X depends only on the spatial coordinates and T only 
on time. What Parker discovered is that this is possible on FRW spacetimes and 
represents a very handy tool for the analysis of the dynamics. A thorough analyt- 
ical investigation of the method in the cosmological context was conducted in j3], 
where an abstract functional analytical eigenfunction expansion was introduced as 



December 12, 2012 2:19 WSPC/INSTRUCTION FILE Bulk 



3 

a methodological background, and precise methods were suggested for the mode 
decomposition of regular solutions on FRW and ultrastatic spacetimes. The theory 
of the method does not seem to have been developed any further until nowadays. 
In particular, the following questions remain open. What are the precise limits of 
applicability of the mode decomposition by means of separation of the time vari- 
able? How many different decompositions are possible for the same geometrical 
setup? When and how can the decomposition be extended to weak (distributional) 
solutions? 

In the first part of this work satisfactory answers will be given to these ques- 
tions and establish a unified framework for the method. Our geometrical setup 
will be a finite dimensional vector bundle over an arbitrary globally hyperbolic 
Lorentzian manifold furnished with a pseudo-Riemannian fiber metric and a fiber 
metric linear connection. The field equation will then be given by the field operator 
D = -|- m*{x) where is the connection d'Alambert operator acting on the 
smooth sections, and m*{x) will be the variable smooth "mass term" (possibly in- 
cluding a coupling to scalar curvature) to which mild assumptions will be imposed. 
This seems to be the most general setup of a (symmetric) hyperbolic linear field 
on a curved spacetime, and covers most practical situations in the cosmological 
context. 

The results can be briefly described as follows. Precise geometrical necessary 
and sufficient conditions are obtained for the mode decomposition of smooth so- 
lutions by time separation to be realizable. This mode decomposition is given, as 
perhaps expected, by the time dependent Fourier transform, and is shown to be ba- 
sically the only such mode decomposition possible. The decomposition is extended 
to all distributional solutions in a natural manner given that there exists a choice 
of modes fulfllling certain regularity conditions. The conditions become fully ex- 
plicit once one has a Paley- Wiener theorem for the spatial Fourier transform, i.e., 
a precise description of the Fourier image of the test functions space. This is the 
situation in FRW spaces. Moreover, it is shown that if the bundle is analytic and 
the dynamics of the geometry sufficiently rigid (precise definitions are given) then 
the conditions are satisfied regardless of the harmonic analysis involved. It turns 
out that that the mode solutions of non scalar fields under certain circumstances 
experience infrared instability periods not known for scalar fields; the author is yet 
not sure about the physical essence of this phenomenon. Apart from this, tradition- 
ally infrared integrability issues arise when integrating modes over the spectrum 
Spec(il'i;) of a Schrodinger operator Dy, if Spec(Ds) includes the eigenvalue 
Here this question of integrability of modes over the spectrum is settled by show- 
ing that the suitably chosen mode solutions remain well under control even at non 
positive spectral values. As an example of application and as a byproduct the ex- 
plicit formula of the propagator of the field in the Fourier space is found, which 
generalizes one obtained in [5]. 

Another advantage of geometric symmetries is the possibility to apply harmonic 
analysis. This is particularly true for the cosmological models where a rather rich 
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group of spatial isometrics is imposed. Then the spatial sections of the spacetime 
can be considered as homogeneous spaces, and the spatial Fourier transform can 
be investigated in much more detail with many explicit consequences. These prop- 
erties then can be dragged to the time dependent Fourier transform and hence 
to the mode decomposition. Of particular interest in the quantum field theory on 
cosmological spacetimes are the spatially symmetric (invariant) states, of which 
the 2-point functions are bi-distributional solutions of the field equation which are 
invariant under the full isomctry group. In [5 the Fourier image of the isotropic ho- 
mogeneous states of the Klein-Gordon field on FRW spacetime has been obtained, 
under an additional continuity requirement which has no clear physical interpreta- 
tion. In contrast the most general form of invariant scalar bi-distributions on M!^ 
has been obtained in [6] using a nice technique. What appears to be missing is 
a generalization of these results to sufficiently many homogeneous spaces so that 
practically all cosmological situations are covered. The harmonic analysis of FRW 
symmetry groups is well known since long, but strictly speaking isotropy is not 
as fundamental in cosmology as homogeneity, and one is also interested in cosmo- 
logical models which are only homogeneous (Bianchi models) or in addition only 
partially isotropic (LRS models) . The isometry groups of these spaces are described 
by Bianchi groups with their quotients and semidirect extensions (in case of LRS 
models). Some of these groups are solvable, others are semisimple, with finite or 
infinite center. Therefore it is not easy to establish a unified harmonic analytical 
approach for all cases, although one has to admit that Kirillov's theory for the solv- 
able groups and Helgason's theory for semisimple groups together would cover the 
majority of situations. To obtain a unified theory one can adopt abstract harmonic 
analysis. This beautiful branch of mathematics allows to obtain many results in 
an admirable generality. However, apart from compact groups, it is not completely 
clear how to relate the abstract group Fourier transform with the eigenfunction 
expansion of the invariant Laplace operator. At least there seems to be no unified 
exposition of these techniques applied in the cosmological context. 

The aim of our second part will be to put together some tools from harmonic 
analysis which are adequate in cosmology, and to obtain useful results with their 
help. First a short survey on homogeneous bundles in general, and on semidirect 
homogeneous bundles in particular is carried out, and it is shown that they cover 
the vast and the most important majority of the realistic cosmological structures. 
Next, the abstract harmonic analytical Fourier transform is introduced with its 
requisites on semidirect homogeneous spaces. Although the abstract Fourier trans- 
form on groups and the representation theory in homogeneous spaces are to be 
considered as well studied and widely known, the abstract Fourier transform on 
homogeneous spaces is not that popular and deserves a better exposition (at least 
we were not able to find a satisfactory one in the literature). Then some properties 
of the Fourier transformed distributions are established. Although some of these 
results may be known to experts in harmonic analysis, we were not able to locate 
them in the for required for our purposes in the literature. And because these prop- 
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erties are needed to obtain our results on invariant distributions they are included 
here with proofs. Next an attempt is made to unify the abstract Fourier transform 
and the usual spatial Fourier transform given by the eigenfunction expansion of 
the Laplace operator. Several remarks are made on this way, which may serve as a 
guideline to completely build the desired correspondence once a particular structure 
is specified. This has been indeed performed for the purely homogeneous cosmolog- 
ical models, which will appear in a subsequent publication. The necessity of such 
an explicit correspondence consists in the ability to transfer the results obtained 
in the abstract setup to the situation with the usual Fourier transform, which is 
far more useful a tool for concrete calculations. Finally, by a generalization of the 
above mentioned technique in [6] the general form of the invariant bi-distributions 
in arbitrary semidirect homogeneous vector bundle is found without any additional 
assumptions on their regularity. It is concluded that the polynomial bound of the 
Fourier transformed homogeneous state as found in [5 is a consequence of the 
imposed regularity requirements. 

2. Mode Decomposition of Hyperbolic Fields 
2.1. Linear hyperbolic fields 

It is generally believed that the forces of nature are described by tensor and spinor 
fields. A geometrical generalization of those are the vector bundle fields, i.e., fields 
as smooth sections of some vector bundles. In general relativity one works mainly on 
a four dimensional Lorentzian smooth manifold (M, g) which is called a spacetime. 
We will be concerned with hyperbolic fields given by a wave equation, hence we 
put an additional constraint on the spacetime (M, g) to be globally hyperbolic, so 
that the Cauchy problem of the wave equation is well-posed. For simplicity only 
linear fields will be discussed here. For the reduction of the Maxwell and Proca 
fields to linear hyperbolic fields the reader is referred to [1],[2]. We summarize the 
basic setup of the the linear hyperbolic fields in the following section. 

Let V be an n-dimensional vector space. Let T ^ Af be a vector bundle with 
standard fiber V and with a pseudo-Riemannian metric {u,v)g. Let further £(T) = 
C°°{T) and 'D{T) = C^{T) be the spaces of smooth sections and of those with 
compact support, correspondingly. Let V be a metric connection on T and 
the associated d'Alambert operator on £{T)- Define the field operator to be the 
normal hyperbolic operator D = +m*{x) acting on f (T), where m* e C°°{M) 
is a generalization of the usual mass term which now can also contain the 
coupling term S.R- Note that because differential operators are support-decreasing, 
D'D{T) C 'D{T). A free linear hyperbolic field (j) E £{T) is a solution of the field 
equation D(j) — 0. 

Being a globally hyperbolic spacetime, M = I x E, where I C M is an interval, 
and for each i e I the hypersurface ~ S is a three dimensional embedded 
Riemannian submanifold, which is spacelike with respect to g and is a Cauchy 
surface in the sense described below. Thanks to [5] one can choose a smooth global 
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time function t and a coordinate atlas such that x — {t, x) — (xq, xi,X2,X3) where 
t £ X and X e S, i.e., Ei are equal t hypersurfaces. The restriction of the bundle 
T to the submanifold will be denoted by 7t = 7r^^(I]t). The spaces of smooth 
sections will be £(7^) = C°°{rt) and 2?{7^) = C^iTt). liit-Tt^Tis the identical 
embedding, then its pullback is the restriction map for objects on T to %■ In 
particular il : £{T) £{Tt) and : 2?(T) 2?(7t) are linear surjective maps. The 
embedding it o o tt~^ : M — > S gives rise to a natural embedding it : TM — TS 
and of all tensor bundles (using the same symbols it, il for different restrictions 
in the spirit of polymorphism should not lead to a confusion). The Riemannian 
metric h on TY, will he h = —il{g), with minus sign here because of the signature 
convention (+, —,—,—). The restriction (V) = = V* is a metric connection 

on Tt- The associated Laplace operator At — A^* is an elliptic operator on £{Tt) 
(so that — At is a positive operator). The restriction of the field operator D to £(Tt) 
will be denoted by D^^ — — At + m*{x). 

An existence and uniqueness theorem [7], [9], [10] for wave operators tells that 
the Cauchy problem is well posed: there exists a bijective linear map 

SiTt) © SiTt) 3 (/o, /i) ^ X/o, /i) e {/ e S{ry. Df = 0} 
such that fo = itif) and /i = iiCVtf), where Vt = V_g_. Furthermore, there exist 

dt 

unique Green's operators : 'D{T) £{T) satisfying DE^ = E^D — idu^-j-^ and 
suppiC^f} C J±(supp{/}) for ah / G V{T). Here J±(A^) with a subset N C M 
denotes the causal future/past of N. Define hy E = E^ — E^ the propagator of 
D, which satisfies DE = ED = 0. Now Sol{T) = ]{£{Tt) © £{T)t) and SoIq{T) = 
j{T>[Tt) © 'D{Tt)) will denote correspondingly the spaces of all smooth solutions, 
and of those satisfying supp{f} n St compact for all t E I, respectively. Then 
E'D{T) C SoloiT). There is a symplectic form on Solo{T): 

a{u,v) = / d^h [{it{u),it{\/tv))g - {i^. {\/ tu) , {v)) g] , Vu,w G SoloiT), Vt G I, 

"'St 

where h = —il{g) is the induced Riemannian metric on St. That this is conserved 
(analogous to a Wronskian in ODE) can be seen by considering the Green's identity 
for u,v & Solo{T) on the regular cylindric region U = (ii;t2) x S C M for any 

ti 7^ t2, 

0= / dHg[{u,Dv)g - {DU,v)g] = 

Ju 

= / d.Hh[{itiu),it{^tv))g ~ {it (ytu), it {v))g],yu,v E SoloiT). 
JdU 

This identity also helps us along with Green's operators to find the explicit form 
of the map j. Given any v £ SolaiT), f G ViT), we apply it two times; once 
for the pair v,u ~ £'+(/) on the region [/+ = (— inf{I};t) and once for the 
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pair v,u ~ E^{f) on the region = (i;inf{I}). Summing up the resulting two 
identities and using the support properties of we arrive at 

/ d^ig{v,f), = cJ{v,E{f)). (2.1) 

JM 

We see that the functional cr{v,E{.)) : I?(T) — > C actually is given by a smooth 
integral kernel, which equals v. Thus we can write symbolically 

j(/o,/i)[y] = / d^ih [{h.VtE{y))g - {fi,E{y))g] , yfoji G V{Tt), yeM.tel. 

For full details of this last computation the reader is referred to illj, where the 
argument is given for 1-forms, but is readily applicable to our more general case. 

Proposition 2.1. The operator E : I?(T) SoIq{T) is surjective. 

Proof. Let v E SoIq, and let Ky — suppw n ([0, 1] x E) be the compact region of 
its support between times and 1. Let further x G £{M) be a smooth function 
which equals 1 for t < and for t > 1. Denote v~ = —vx and v'^ — v{\ — x), 
then V = — . Let — Dv^ , then supp/^, C is compact, hence fy G 'D{T). 
The equation f^, = Dv^ has a unique solution with past compact support, and it 
is given by w+ — E^fy. Now Dv^ = —Dv + Dv^ — /„, and similarly — E~ fy. 
Then v — E^ fy — E^ fy — Efy. The arbitrariness of x reflects the non-injectivity 
of □ 



2.2. Spectral mode decomposition 

Henceforth we will use nomenclature introduced in the Appendix A without special 
notice. Consider the operators D : V{T) ->■ V{J) and D^, : V{Tt) V{Tt). If 
m*{x) G M everywhere on M, then by the virtue of Green's identity D and Ds^ 
are formally self-adjoint with respect to the inner products (, )m and (, We will 
not need the self-adjointness of D in the current work. The constructions below 
will pertain mainly to i'st- The conditions on m*[x) for D-^^ to have a self-adjoint 
extension can be found in |12j . We moreover require that the operator Ds^ be lower 
semi-bounded. In practice will be mainly interested in cosmological models, where 
m*(x) = m*{t) is a function of time only, so that no problems arise. Below we 
assume self-adjoint extensions for both D and D^^ , but for D this is only symbolic 
and targets simply at coherent notations. 

Let D and be extended to self-adjoint operators on L^{T) and L^iTt) re- 
spectively. In the rigged Hilbert spaces i],[13],[14] V{T) C L^{T) C V{T)' and 
I?(7t) C L'^{Tt) C ^{Tty operators D and £>St possess complete systems of eigen- 
functions {up} and {Ca} satisfying 

Dup = pUp, Up e ^{Ty, p e M, 
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Denote by ^{Typ and 'D{Tt)'x linear spaces of eigenfunctions corresponding to 
p and A, respectively. Furthermore, there exists an isomorphism 

L^(Tt)^ r dv{\)H{\), (2.2) 
Js. 

where 

-DsJ_fi-(A) = A, 

and dv{\) is a positive measure. Each H{X) is continuously embedded in 'D{Ttyx- 
The eigenfunction expansion of D will be the map 

V{T)3f^fpe{V{rypy*,peR, 

{X'* denotes the space of continuous antilinear functionals on the space X) where 
fp is defined by 

fp{up) = upif), V/ e V{T), Up e Viryp. 

(Here we defer a little from Gelfand's notations who puts fp{up) — Up(/).)The 
expansion of on is constructed similarly. Note that Z^s^ is an elliptic 

operator, hence ^^{Ttyx C £{%)■ 

If each 'D{Ttyx finite dimensional (eigenvalue A has a finite multiplicity N\), 
then 

H{X)=V{Ttyx, dimH{X)^Nx. 

Choose {C\.i}il^i be a an orthonormal basis in 'D{Ttyx (orthonormality understood 
in H{\)). Then {V{Tt)'xy* ^ by the bijective linear map 

{/. = /(CA.)}t\,V/e(I?(7;)',)'*, 



/(Ca)-/(^E^.Ca,»^ 



where each fi e C. In particular, if fx is the mode expansion of / G ^{Tt), then 
the map 

ViTt) 3 f^fx^ {fxA e r di^{X)C^' (2.3) 

Jm 

will serve as a Fourier transform on 'D{Tt)- Define 

Spec{Dj:^\ = suppjdi/}, 

and 

S = {(A,i): A e 5pec{Z?sJ, « = 1,...,^a}. 
Define the spectral measure on S as 

dp{a) — dv{\) X d^{i), 
where djj is the counting measure. The map fEa l2.3p can be reformulated as 
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Then the formula fEa l2.2[) arises a Plancherel formula 

(/, / d^i{a)f{a)g{a), 
and a Peter- Weyl (or Fourier inversion) formula 

f[x) = f dii{a)f{a)Ux), (2.4) 

which holds in the L^-sense so far. In our cases of interest this convergence will be 
in the compact topology. 

However, if 'D{Tt)'\ is infinite dimensional, more delicate tools are needed to 
obtain a Fourier transform with desired properties. Such tools naturally include an 
investigation of symmetries of the underlying geometrical structure, and this is the 
subject of the harmonic analysis. We will often use the formal structure fEq |2.4p 
without mentioning a concrete realization, assuming that this is possible. For the 
cases of our interest we will indeed find a realization by means of adapted Fourier 
transform. 

In the theory of Fourier transform, and in particular in the Euclidean case, the 
Paley- Wiener theorems describe the functional analytical structure of the image 
F\D[Tt)\ of the test function space under the action of the Fourier transform. This 
description is very useful when analyzing the situation in the Fourier space. Unfor- 
tunately there is no (at least known to us) general Paley- Wiener argument valid for 
any Fourier transform arisen in this manner, and the proofs of the existing ones are 
rather structure-specific. In applications we would like, however, to obtain results 
which arc valid in a large variety of cases, and therefore we will introduce a no- 
tion of 'conventional' Fourier transform which consists of a number of assumptions 
pertaining to the analytical properties of a given Fourier transform. Some of our 
later results will be valid under the assumption that the eigenvalue expansion of the 
operator At has at least some of the properties of a conventional Fourier transform. 
One says that a good definition is an assumption of a theorem. In this sense the 
following is not a good definition as we will not manage to use all properties in this 
work. However it seems feasible that these properties will become useful for several 
applications in quantum field theory. 

Definition 2.1. A Fourier transform J- given by the eigenfunction expansion 
against a complete system {Calags will be called conventional if 

(i) The Fourier space (or momentum space) E is a manifold consisting of n = dim V 
components, E = IJ"^]^ S*, and each component is either a discreet set or an 
(not necessarily connected) analytical manifold 

(ii) The eigenvalue A (a) is an analytic function on E 

(iii) The range F\V{Tt)\ is a subspace of the space of analytic functions f{a) on E 
with rapid decay in A 

(iv) There is an involution a — > —a on E such that C-a = Cq- 
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Note that it follows A(— a) = X{a). In later sections we will give harmonic 
analytical justifications for such a 'conjecture'. This conjecture is anticipated, in 
particular, for all cosmological models. Moreover, being true for FRW spaces, it 
can be proven also for Bianchi I- VII spaces (to appear in a future publication). 

Further we will be mainly interested in the space of weak solutions of the field 
equation, I?(T)o, and will try to find a convenient characterization of it. In partic- 
ular we will be looking for a complete system of solutions {uq} spanning I?(T)o 
and being in addition well handled (i.e., smooth, explicit etc.). One means of doing 
this is to look at a subspace of 2?(T)o which consist of solutions f(x) = a{t)b{x), 
a e C°°{I), b £ £{%)■ Then under fortunate circumstances the field equation breaks 
apart into two lower dimensional elliptical eigenproblems, which are much easier to 
deal with. Which are those circumstances and whether such solutions span 2?(T)o, 
and related questions, are the matter of the problem of variable separation. In the 
next sections we will find out in which cases this is possible and how to perform it. 

2.3. Separation of variables 

As discussed above, we would like to span the space 2?(T)o of weak solutions of 
the field equation by a family of easily computable smooth solutions {ua}- In this 
section we will see when and how one can perform this for the smooth solutions 
Solo{T). The necessary requisites for this will be predominantly geometric require- 
ments. In the next section we will show that under additional functional analytical 
assumptions the procedure can be extended to 2?(T)g in a natural way. 

Definition 2.2. Let be a topological (complex) vector space with closure S ^ S, 
9Jt a measure space with measure dm. An 9Jl-measurable family {ualagOT of linearly 
independent elements Ua € S will be called a complete or spanning system for S if for 
any v € S there exists a unique (modulo null-supported functions) 9Jt-measurable 
function aV ^ R (a" : OJl ^ C) such that 



with the integral converging in this topology. 

Note that uniqueness is implied already by the linear independence. We will 
always take 9Jl to be minimal, i.e., there exists no subset A C VJl with m(j4) > 
such that a"" {A) = for ah u e S". 

£{T) is a closed topological vector space with the topology of compact con- 
vergence, and Sol{T) and SoIq{T) arc linear subspaces. 'D{Ty is a another closed 
topological vector space with its distributional topology, and 'D^TYq is a linear sub- 
space. Through the natural embedding, £{T) is also a subspace of I?(T)', and the 
compact topology of £{T) is stronger than the distributional topology of I?(T)' 
(|TS], Chapter XVII). Therefore if {ua}a&m is a complete system for Sol{T) with 
compact topology, then it is such also with the weaker distributional topology. 



V = 




'Q 
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A spanning system {ua}aem of Solo{T) of the form Ua = T^Xa, where G 
C°°{I) and Xa G such that Dua = 0, will be called a complete (time-)variable 

separated system of solutions (or shorter, separating system). 

We will assume that a Fourier transform on T>{Tt) is specified by means of 
the spectral decomposition of as described in the previous section. The system 
of eigcnfunctions {C^} of A^, with the Fourier space St and the spectral measure 
d^{a) on it, provide a spanning system for T>{Ti) by means of the Fourier inversion 
(or Peter- Weyl) formula. Below we will come across the question of a spectral 
theory of formally non-self-adjoint, i.e., asymmetric differential operators of type 
a{x)D^^ . As a rule, the cigcnfunction problems of asymmetric (aside from unitary) 
operators are ill-posed, and eigenfunctions do not comprise a complete system, but 
there are rare exceptions. At this point we have to admit the non-exhaustiveness 
of our treatment, as we do not analyze this possibility. We will loosely rule out the 
possibility of such operators to have a well-posed eigenfunction problem. 

A small remark will be useful later in the section. 

Remark 2.1. If {TaXa}aeon is a separating system for Solo{T) with compact 
topology, then for each t Gl, both the families {Ta{t)Xa}aeTi and {Ta{t)Xa}aeM 
are spanning systems for 'D{Tt)- In particular, for each x & Y,t, the family 
{Xa{x)}a£m contains a (possibly redundant) basis of V. 

The assertions are relatively obvious in the view of the fact, that the restriction 
maps i( , i( o V( : Solo{T) — )■ 'D{Tt) are surjective, and hence a spanning system for 
Solo{T) must give a spanning system for the Cauchy data V{Tt) (B 'D{%) on S^. 

Proposition 2.2. Let {TaXa} be a separating system for Solo{T). Then for dxa- 
almost each a € 9H there exists at least one ^ 3 (3 ^ a such that Xa = X^. It 
follows that Ta and are linearly independent. 

Proof. First we note that whenever X^ = X^ for a 7^ /3 it follows that and 

Tfj are linearly independent, otherwise the linear independence of TaXa and T^Xp 
would be violated. This proves the second statement. For any u S Solo{T) we write 



For a pair 77 = (?7o,?7i) G (C^) set (I'o.w?) = (%uo,'7iWi). Then = j{vl!,,v'^) 
will be a solution v'^ G Solo{T). 

Let SOti , 9712 C be the subsets for which Xa = Xp happens for no and 
for one or more pairs a,/3, respectively. Because {TaXa} is measurable, the sets 




The Cauchy data {uo,ui) of u on St are given by 
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mi,M2 cM are measurable, and 97li +0^2 = Fix a time tel. By Remark 
the fields Ta (t) and Tq, (t) vanish dm-almost nowhere on 931. For a E DJli , a ^ /3 G dJl 
we have Ta{t)Xa and T/3(t)X;3 are linearly independent in S{Tt)- For any a G 2)t2 
there is at least one /3 G 93t such that P ^ a and Ta(t)XQ and Tfi{t)Xp are linearly 
dependent, and all such (3 lie in 9Jl2- It follows in particular, that whenever 

dm{a)a''{a)Ta{t)Xa = 

m 
then 



dm{a)a^{a)Ta{t)Xa= / dm(a)a"(a)r„(t)X„ = 0. 

The same is true also with Tq, (t) instead of Ta (t) 

Let u € SoIq{'T) be a nonzero solution, and choose time t such that ui ^ on 
Ej. Such a choice is always possible because static functions do not solve the field 
equation. Set rj = (0, 1), then will be a nonzero solution with Vq = on St. Then 



0= / dm(a)a" (a)T„(t)X„ = / dm(a)a" (a)T„(t)X„, 

and because Ta{t)Xa are linearly independent for a e 93ti, it follows a^'\a) = 
on 93ti . But on the other hand 

0=1 dm{a) la"" (a) - a"(a)j ta{t)Xc, = [ dm{a) U"''(a) - a"(a)j fa{t)Xc, 



mil 



from where it follows a"'' (a)— a" (a) = and hence a" (a) = dm-almost everywhere 
on 93Ti. But u was arbitrary, and because DJl is minimal it follows m(93li) =0. □ 

Remark 2.2. Let two equations f{t) + F{t)f{t) + G{t)T{t) = and f{t) + 
H{t)T{t) + J{t)T{t) = have two common linearly independent solutions T{t) and 
S{t). Then by Liouville formula the Wronski determinant det W[T, S]{t) evolves by 

detW^[T,5](i) = dctW[T,S]{0)e-J"'^^^<-^'^ = detW[T, S]{0)e- Jo 

hence F = H and thereby also G = J. 

Proposition 2.3. The solution space Solo{T) admits a separating system if and 
only if there exists a covering of T by local trivializations such that the following 
local conditions are satisfied everywhere (metric g is time-separated): 

fi) 9oa — 300(^)7 the metric component goo depends only on time 
(ii) the expression Y^^i j=i 9^'' {x)^^ {x) is a function of time only 
(Hi) the connection 1-form T and Christoffel symbols F*^ satisfy 
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3 



k=0 



Tq = Toit) is a function of time only 

(iv) the eigenfunction problem of Dy,^ on different St can be adjusted, so that all 
Tit cLre isomorphic and the eigenf unctions Ca — Ca o-re time-independent. 

Proof. Throughout the section we will work exclusively locally, i.e., in a local 
trivialization tt~^(U) U x V , U C M. Thus we identify the sections in a bundle 
having a typical fiber with functions in C°°{U; 5"). We will not keep the flag U in 
this section but will always understand objects as restricted to U. 

The d'Alambert operator on £{T) has the following local expression in terms 
of the connection form coefficients Ti and Christoffel symbols F*^^- , 

92 



3 

E 

i,j=0 



9- 



fc=0 



' dx'^ dxi 



fe=o 



C2.5) 



and the field operator D locally looks like 

3 02 



^=E 



52 



9- 



i,j=0 



3 

4=0 



where A*,_B € C°° {U, End{V)). To achieve a time separation we need to choose 
a coordinate atlas such that everywhere = for i > 0. Then the operator D 
locally breaks apart into two differential operators, D = Dt + D-^^, where 



r.00 



dt^ 



at 



and 



D^t= y . 9- 



3 

E 



d d 

dx^ dx^ 



E^^#- 

^ dx' 

i—1 



I* = -A, 



is the restricted field operator defined earlier. B^^B^ e {U , End{V)) are to be 
seen explicitly from (Ea l2.5p . 

^Necessity: Let {TaXa} be the separating system system. Then 

DTo,{t)X^{x) = [Dt + i?sjr„(t)X„(f) = f,(t).9™(a;)X„(f) + 

+r„(t)A°(a;)X„(f) + T,(<) + D^^] X^ix) = 0. (2.6) 

That the metric signature is definite it follows that never vanishes. We find 

a family of second order linear homogeneous differential equations 

f„(i)5°"(x)X;(x) + T„(t) {A%x)X^{x)y + T^{t) {[B%x) + D^,] = 

parameterized by the spatial coordinates G S and fiber indices i = l,...,n. By 
Proposition 12.21 we know that there exists a f3 ^ a with Xjj ~ Xa and , 
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linearly independent. This means that all these equations share at least two linearly 
independent solutions and Tg. If for some x and i, X'l^{x) = 0, then the existence 
of two linearly independent solutions for the resulting first order equation means 
that 

{A%x)Xa.ix)y = {[B%x)+Dj:,] X^ix)y - 0. 

Otherwise, by Remark 12.21 we find that there exist functions Fa, Ga G C°°(I) such 
that 

{A%x)X^{x)y ^ g'"{x)Fa{t)Xl^{x), {[B%x) + D^,] = g°"(x)G„(OX;(f). 

In both cases we establish that 

gooix)A''{x)Xa{x) = Fa{t)Xa{x) (2.7) 

and 

fMx) [S"(a;) + X^{x) = G„(t)X„(f). (2.8) 

Thus for each t Cz T, Xa-s must be nothing else but the joint eigenfunctions of 
the operators gooix)A^{x) and goo{x) [B°{x) + -DeJ corresponding to eigenvalues 
Fa{t) and Ga{t), respectively. The operator goo{x)A'^{x) is simply a matrix, and at 
each point x G M has at most n independent eigenvectors. By Remark |2.1[ Xa{x)- 
s span V, and thereby {Xa}aGm contains bases of all eigenspaces of goo{x)A^{x). 
From (Eq l2.5p we find 

3 

goo^o = 2ro-goo 51 ff'-T"^' (2-9) 

and 

^3 3 

fc=l i,j=0 

Now turn to the eigenfunction problem (Eq l2.8p . As discussed above, for 
this problem to be well-posed it is necessary that the differential operator 
300(3^) [B^{x) +£>St] is at least formally self-adjoint. But this is possible only if 
goo{x) — goo(t), thus we have obtained the condition (i). Let us switch to an atlas, 
where the time function t is redefined such that 500 (t) — 1 (this step is not crucial, 
but only for convenience). It follows, that 

rgo = 0, Vfc > 0, 

so we obtain 

3 

AO = 2ro- 5] g*^r°^-, (2.10) 

B" = ^ro + rg- ^5'^r° ro. (2.11) 
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Combining fEa l2.7|) and fEa l2.10() we see that {Xq}-s are the eigenvectors of Tq, 
and these eigenvectors are independent of t. Hence they are also the eigenvectors 
of §1^0, and thus by fEa l2.11( ) and are simuhaneously triangularizable, 

for some i7„ e C'°°{M). We note that 



^ 1%. 
2 dt 



and denote 



P{x) = - ^ 5^^(x)r°.(x) = i ^ 



Now fEa l^TT)) and (EqU^ tell us, that for each i e I the operators A° and + 
have a common system of eigenfunctions spanning 'D{Tt), and therefore must 
commute. 



This requires 



and 



[A", D^^ +B"]u= [AO, D^^] u = 0,Wue V{Tt). 
AO(x) =2ro(x)+P(x) = 

3 

^g'-'"[ro,r,] = o,Vj>o, 



ro,ff + r.r,-;^r^r. 



k=0 



= 0, 



exactly as the statement. Similarly, that operators B^ and 5° + have the same 
eigenfunctions implies, that [i?o,Z?sJ = 0, which on its turn requires B'^ix) = 
B°{t), and thereby P{x) = P{t) and ro(a;) = ro(<). Thus we have proven parts 
(ii) and (iii) of the statement. It follows further, that Ha{x) = Ha{t), and thus the 
eigenfunction problem (Eq l2.8p becomes 

D^^X^ix) = (G„(t) - i/„(i))Xa(x). 

This is exactly the eigenfunction problem of , whence we conclude, that neces- 
sarily 



{Xa}a€m C {({} 



Therefore 
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where 

A„(t) = {AeM:X„e2?(7^)U- 

Now (Ea l2.6p becomes 

f^{t) + F^{t)f^{t) + G^{t)T^{t) = 0, (2.12) 

which is the mode equation for the mode T^. We have two spanning systems for 
^{Tt): {Xa}aem and {CalagSji and hence in each eigenspace ^{Ttyx we can choose 
a basis from {Xa}aew- Thus a complete eigenfunction system can be chosen among 
{Xa}ai£9yti proving the (iv) statement of the proposition. We are complete with the 
necessity. 

■^Sufficiency: Suppose all the points of the statement are satisfied. Then, as we 
have seen above, by (iii) A° and are functions of t having the same eigenvectors, 
and moreover, commute with A^. It follows that the actions of A'^ and 5° preserve 
'D{Tt)'\, and thus by a Gramm-Schmidt operation the representatives Ca can be 
chosen such that they are eigenfunctions of and . Thus each Yj\, and thereby 
the entire S, decomposes into n components corresponding to the eigendirections 
of AO, 

n 
i=l 

For spatially homogeneous spacetimes discussed in later sections we will give a more 
conceptual justification of such a subdivision in terms of the representation theory. 

Let for each a € E choose a mode solution of fEa l2.12( ) arbitrarily (strictly 
speaking, not completely arbitrarily, but such that and Tq, are linearly indepen- 
dent) and consider the union of two systems 

Choose any € Solo{T)- Then for each t e X the restriction i'l{(f>)[x] — cl){t,x) G 
'D{7t) can be Fourier expanded as 

(l){t,x) = J dn{a)${a;t)<:a{x) (2.13) 

with the integral converging in L^(E,/i). Hence we can differentiate under the in- 
tegral, 

D(p{t,x) ^ / dii{a)D (j){a;t)(a{x) = dfi{a) (j){a;t)+ 



s 



E 



+F„(i)(^(a;t) +G„(i)0(a;t)J Ca(^) = 0, 

where for convenience we again reparameterized t to get 500 = 1- Thus 4>{a;t) is a 
solution of the mode equation. All solutions of the ordinary second order equation 
(Eq |2.12]) are smooth and comprise a two complex dimensional space. 
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Inserting this into (Eg 12. 13]) we finally arrive at 

which exactly means, that is a separating system for SoIq{T). (For com- 

patibility with the definition one can concatenate Ua and Va to a single function on 



The assertion of this proposition can be interpreted as follows. If a mode de- 
composition in a reasonable fashion exists for SoIq(T') then it is basically the mode 
decomposition given by the time dependent Fourier transform which we will define 
a few paragraphs later. 

As a supplement to the proposition we make a few remarks. Let gij = g{di, dj) 
and hij = h{di,dj) be the matrices of the metrics g and h, correspondingly, in a 
local chart, and <Tk{h) the eigenvalues of the symmetric matrix h. 

Remark 2.3. The condition (ii) of Proposition 12.31 is equivalent to 

detg{x) = -gooit)deth{x) = -gooit)'Jiihy2ih)c73{h) = -gooit)e^ ^« ^^''Uethoix), 

where detho{x) € C°°{Y.) is a positive smooth function (the notation will become 
clear later). 

The assertion follows from the combination of condition (ii) with the Laplace's 
formula, 



Because V is a metric connection, the restrictions of the previous proposition 
imply restrictions on the fiber metric (, )g. In case of a tensor bundle of rank (m, n) 
with Levi-Civita connection, coefficients F^ are expressed in Christoffel symbols and 
the fiber metric is induced by the spacetime metric, thus the restrictions fall onto 
the spacetime (Af , g) . 

Corollary 2.1. Let a local moving frame be chosen, such that the metric is 
represented by the matrix g. Conditions (Hi) of Proposition [273\ imply the following 
restrictions on g.- 



bundle of rank (m, n) to allow for seprataion it is necessary that the spacetime 
metric be represented by a matrix 



the disjoint union E U S.) 



□ 



-det.g = det5-Tr[r'7^5]- 





where ho and B are matrix valued smooth functions. 
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Proof. Locally the conservation of the metric V(,)g = can be written as 

where Ti are the matrices of the connection coefficients in the chosen frame. In 
particular, for i = we have 

where Fq = ro('^) was used. The solutions of this equation are of the form 

g{x) = %^m°{x)^{t), 

where 

^{t) ^ e^o dtToit') ^ (2.14) 

and g'^(a;) is a smooth symmetric matrix field on S. 

Now if we identify the tensor space (TpM)™ with an 4"+"' dimensional vector 
space V using a suitable bases, then each matrix F^ will be a 4"+™-! x 4"+™-i 
matrix of blocks, with blocks being the Christoffel symbols for contravariant 
indices and — F^ for covariant indices. Fq = Fo(t) means Fq = Fo(i). With our 
time-separated metric we have 

g=l®-h. 

One can find 




for some smooth 3x3 matrix A{t). The solution is 

h{x) = ho{x)e''^odt'Ait') ^ ho{x)B{t), 
for smooth symmetric commuting matrix fields ho{x) and B{t). □ 



Now the notation det ho of Remark 12.31 becomes clear, and we see that 

detB(i) =e2/o'imt') 

for a tensor bundle. Note that for the scalar field conditions (iii) are trivially satisfied 
and do not restrict the spacetime. 

Remark 2.4. For the volume form measure dfXfi on Ei we have locally 
d^hix) = \J det h{t,x)dx^dx^dx^ . 

By Remark 12.31 we have 

det h{t, x) = fo det ho{x), 
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hence 
where 

dfihoi^) = detho{x)dx^dx^dx^ . 

Henceforth by stating that a mode decomposition of SoIq{T) exists we will mean 
that the assumptions of the Proposition 12.31 are satisfied and the corresponding 
covering is chosen. We are ready to formulate precisely the time dependent Fourier 
transform. Note that although C,a are i-independent, the spatial metric h and the 
fiber metric (,)g depend on t, and Ca are not orthonormal with respect to the 
measure d^,h for all t simultaneously. At this point we appoint once and forever to 
normalize Ca such that they are orthonormal at t = 0. Or equivalently, they are 
orthonormal with respect to the measure d^.^^ of Remark 12.41 and the fiber metric 
gO of Corollary [HH 

Definition 2.3. For / G 2?(T) we define the time dependent Fourier transform 

= F[f(t, .)]{a) e by 

-^[/(^) •)](") = / dfIho{x){Caix)J{t,x))gO. 

Here we note another important corollary, which will be useful later. It will give 
the time dependent Plancherel formula. 

Corollary 2.2. Suppose the assumptions of Proposition \273[ are satisfied, and the 
corresponding covering is chosen. Then for all f € T^{T) 

(Ca,/(i,.))s, .)](«) 

and the time dependent Plancherel formula for the time- dependent Fourier trans- 
form is given by 

= ^d^iia)Ut)T[f{t;J]{a)T[h{t,.)]{a), 

where 

Proof. By Remark [2:41 

(/(t, .)Mt, 0)2, - / d^,h{s){f^t, .),h{t, .)), = 

At the same time by Corollary [O] we have 

(/(t, .)Mt. ■)), = m)fit, .)Mmt, •))bo- (2.15) 
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Because we have normalized Cq with respect to dnho and 0°, the conventional 
Plancherel formula holds for them, 



d^ihMifit,■),Ht,.))^o ^ / d^^{a)T[f{t,.)]ia)T[h{t,.)]{a). 
St "'s 



Combining these three formulas we find 



is 

Meanwhile 

iCc Jit,. )h,^e^o Pin f dfiHo{x)mt)(:c.,ntMt,.)),o. 



St 



By definition 
whence 

{C»J{t, .))st = fo ^*'[^»(*')+^(0]^[^(i)/(i, .)](„). 

Finally 

= / d^,{amf{t,.)]{a)m{t)U^) = / dfi{a)^[f{t,.)]{a)eifo^''i^-(t')-P(t'^^Ca, 
Js Ji; 

thus 

The assertions now easily follow. □ 

At last we compute the spectra of operators and for the tensor bundle to 
find the functions and Ha- In view of Corollarv l2.1l thc function P{t) becomes 



Then 

Spec{A^} = 2Spec{rQ} + P{t), 

Spec{B°} = {cT(t) + a^(t) + a{t)P{t): a e Spec{ro}} . 

As a useful example we calculate these spectra for the scalar and 1-form fields on 
miiformly expanding (e.g., FRW) manifolds, 

ds^ = df - a^(t)da'^{x). 

Here the matrix B(t) — a?{t)l, and hence 

A{t) = ^lna(t)l = i/(t)l. 
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and 

For scalar case n = m = and we have 

SpeciTo} = {0}, 

thus 

Spec{A°} = {3H{t)}, Spec{B°} = {0}, 
as well known. For the l-form case, m — 0, n — 1, we have 

Spec{To} = {0,-Hit)}, 

and thereby 

Spec{A^} = {3H{t),H{t)}, Spec{B°} = {0,-H{t) - 2H^{t)}, 

where the first members are similar to the scalar case and represent the scalar 
modes, but second ones represent the transversal and longitudinal modes. 

As we have seen, for the separation it is necessary that the evolution of the 
metric be represented by linear transformations. If the connection also satisfies such 
a condition in a suitable sense, than the operator D^t is essentially the same at every 
t up to some scale factors. (Maybe the condition (iii) of the main proposition already 
implies such a restriction on the connection, but we are not sure yet.) This will be 
the case for all our bundles of interest, and it will provide analytical advantages. 
To summarize what we expect precisely we give the following definitions. 

Definition 2.4. We will say that the operator Dj:^ has a strictly uniform 

spectrum over time if there exists a lower semi-bounded function Lu{a.) on E, 
a positive smooth function C{t) > and a smooth function m*{t) such that 
Aa(i) = oj{a)C{t) +m*{t), or equivalently, the expression 

^ln|A„(0-m*(i)| 

does not depend on a. 

This is a rather strong condition. It basically requires that the eigenspaces of 
DY,t coincide for different t up to an overall shift, and that eigenvalues be linearly 

proportional. Such a property would be very comfortable, but it does not hold for 
some models of our interest. In particular, it does not hold for the Bianchi I model 
with distortions. Hence we will derive some of our results under a milder restriction 
which holds at least in all cosmological situations where the spectral theory is 
explicit enough so that it can be checked. (The explicit spectral theory of the scalar 
field on Bianchi I- VII spacetimcs will appear in a subsequent publication.) 

Definition 2.5. We will say that the operator Dy.^ has a loosely uniform spec- 
trum over time if 

d 



^^\n\K{t)-m*{t)\ 



< Cn, ytGn,aGi:, 
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for any compact interval TZ dX, and for some < Ctj G M and a smooth function 
rri*{t). 

If the Fourier transform is conventional, then it will be natural to require that 
w be an analytic function on E. 

2.4. Some properties of the mode solutions 

In this section we investigate the equation (Eg 12. 12]) and obtain some useful prop- 
erties of the mode solutions Tq . The mode equation is 



and \a{t) are the eigenvalues of the operator D^^ = — At + m*{x). Note that Ga 
may become null or negative for some rates of expansion. This corresponds to the so- 
called positive back-reaction in a linear system and results in exponential solutions. 
This is an interesting phenomenon appearing in non scalar fields (for scalar fields 
Ha = 0), and its significance is not yet completely clear to us. To understand it 
one could, for instance, track its influence on the energy-momentum tensor etc. 
It is not obvious that this is really a physical infrared instability, because it may 
occur for the co-vector field but not for the vector counterpart, for instance. It is 
also worth mentioning, that for the co- vector (1-form) field the introduction of a 
conformal coupling precisely cancels this instability. It seems plausible that for each 
field there is a choice of the coupling constant which compensates this bad infrared 
behavior. We say infrared, because Xa[t) attains arbitrarily large positive values at 
any t, thus on an unbounded subbundle of i? x E, Ga is positive. 

We will make this more explicit under the assumption, that has a strictly 
uniform spectrum. Then the function C{t) is uniformly bounded from below, and 
the functions Ha and to* are uniformly bounded from above on any compact in- 
terval TZ. On the other hand a; — ^ -f cx), hence is suffices to choose w large enough 
to make Ga = Ha + Cuj + to* > 0. 

Fix a component E' and write H — Ha, F = Fa and I = la for all a S E'. 
Define a new variable 



which is in a smooth monotone bijective correspondence with t. The inverse function 
will be denoted by t(s). Regarding all the acting functions of t as functions of s we 
obtain 



Ta{t) + Fa{t)Ta{t) + Ga{t)Ta{t) = 0, 



where 



Ga{t)^Ha{t)+\a{t), 




Ta{s)+Ka{s)Ta{s)^Q, 



(2.16) 



where 



Ka{s)^\Ga{t)^^o'^^nr) 



. t=t{s) 



Ga{s)l\s). 
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This is a time dependent harmonic oscihator equation, to which the results in the 
appendix apply. 

Remark 2.5. Note that the Wronski determinant of two solutions Q, R 

det W[Q, R]{s) = Q{s)R{s) - Q{s)R{s) = const 

in variable t becomes 
dt 



det 



W[Q,R]{t) ^ — [Q{t)R{t) ~ Q{t)R{t)'^ = I{t) {Q{t)R{t) - Q{t)R{t) 



— const. 



Applying Corollary [Appendix B.l| to fEa l2.16() for different a we find estimates 
which in principle depend on a in a complicated way. But under the assumption of 
loose uniformity on we will be able to invoke more comfortable expressions. 

Proposition 2.4. Suppose D^^ has a loosely uniform spectrum over time. Then 
for a family of arbitrary solutions T„ of fEa \2.l6\) the following estimate holds 



\Tais)\ < Rn\Ta.m 



S-jz 



max{l, ^UTz+TTzXa{Q)} 
with < Rtz, Stz, T-jz G M and U-jz G M, for any compact interval TZ. 

Proof. Fix a compact interval TZ and for each a G apply Corollary 
[Appendix B.l with Aq(s) — P{s)H{s) + P{s)\a{s). Because Aq, is real, we get 
A-iz{a) = 0. As Xa{s) is lower semi-bounded we have 

Ptz = inf inf Aq > — oo. 

Denote m-jz = mf-jzlPH} and n-jz — inf-jzlP} > 0. Then 

CTzia) > m-jz + n-jz inf Xa > mn + nnpn- 



It follows that K{a) < y^l + \m-iz + n-izPTz\ and 67^,(0;) > 1 + max{0,m7^ + 
n-jzinfTzXa}. Denote M-jz = supT^{\PH\} > and N-jz = supT^{P} > 0. We 
find next 



DTzia) < 1 + IniTz + n-jzPv] + M-jz + N-jz 
Now we observe that by loose uniformity 
|A„(s)-m*(s)| 



sup Xa 
TZ 



In- 



|A„(s')-m^(s')| 



dads In |Aq(ct) — m*{a)\ 



< ITZWNtzCtz, 



hence 



sup |A„ - m*\ < |A„(0) - 7^1-^(0)161^1^^^ 

TZ 

inf |A„ - m*| > \XaiO) - m*(0)|e-l^l^^'^^. 



(2.17) 
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Note that whenever Aa(0) — to*(0) > then it follows by continuity that Aq(s) — 
m*(s) > for all s eTZ. Denote 

A„„ = TO*(0) - niin{0,inf m* • el^l^^-^} - min{0, H^e^'^^^^^}. 

Then from Aq(0) > Xmin it follows Act(O) — m*(0) > 0, m-jz + n-jz inf^ Aq > and 

inf A„ > inf + (A^^fO) - m*(0))e-l^l^^^^ > 0. 
n n 



Now we have that 

en{a) > 1 + x[Aa(0) > A„„;„] (rrm + rm{Mm* + (A,(0) - m*(0))e-l'^l^^^)) 



where the characteristic function x plays here the role of the condition checking. 
From (Eq i2T7l) we find 



supIAqI < supm* + |Aq(0) - TO*(0)| 



n 



n 



whence 



Dn{a) <l + \mn + nnPn\ + Mk + Nn{supm* + |A„(0) - 7^1-^(0)161^1^^^). 

Thus we establish that for Aq, (0) < Xmin 

< 1 + Im^ + nnPnl + Mn + Nnisn^rh* + (|A„„„| + |7?i*(0)|)el^l-^^-), 

and for Aq(0) > X^in 

Dnia) ^ 2\TZ\-/Wi^c^ jy ( J_ + 1 + \mn + nnPn] + + jV^I sup,^ m*\ 



Finally 



— ln(K2+A„) = ^ 
as 



^{P{H + m*)) + ±{P){X^ - 77V) + I\X^ - m*)±\n \X^ - m*\ 
K? + PH + J2Aa 



Denote P-ji = supT^ \ds{P{H — fri*))\ > and Qn = sup^ \ds{P)\ > 0. Again using 
the loose uniformity, for Aq — 771* < 1 



and else 



ds 



ds 



IniK^ + Ac) 



ln(K2 + A„) 



<PTZ + QTZ + iNn)^Cn, 



< 
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Summarizing this all we find that by Corollary [Appendix B.l there exist numbers 



< Rn , S-ji , T-ji e M and U-jz G M such that for a family of arbitrary solutions 
we have 

|T„(.)| < Rn\T^m + J\ rnn '^"^"^'' 

maxjl, ^yu^z + TkAq(O)} 

what was to be proven. □ 

The result can be strengthened under additional assumptions. These are perhaps 
too restrictive, but they appear to be sufficient for some important applications. 

Let Ha = {z e C : < a}. 

Proposition 2.5. Suppose the bundle T is analytic, so that all functions figuring in 
fEa \2.16\) are real analytic functions of s. Suppose further that D^^ has a strictly 
uniform spectrum. Choose the initial data to be Tq(0) = p{u}{a)) and Ta{0) — 
q{uj{a)), where p{i-^), <z(w) are holomorphic functions on Hq for some a > 0. Then 
for each s, Ta{s) — rs{uj{a)), where rs{uj) is holomorphic in uj on Ma and real 
analytic in s, and for any compact interval TZ it holds 

ksHI < Rtz\pH\ + — —-^^===-\q{uj)\, vs e n, 

maxjl, VU-jz + InJi^l 
with < Rn, Sn, e M and Un e M. 

Proof. By strict uniformity we have Xait) = u){a)C{t) + fn*{t), and if the initial 
data depend only on w, then the solutions will also be such. Therefore for conve- 
nience we write 

2L(s) + I^is)iH{s) + ujCis) + m*is))T^{s) = (2.18) 

with Ti^{0) ~ p{i-^) and Ti^{0) = q{io). From the theory of power series it is clear 
that any real analytic function on s(X) can be extended to a holomorphic func- 
tion in some open neighborhood d{s(I)) of s(I). Consider (Eq l2.18| ) as a complex 
differential equation, then for any uj G Ha, by Satz 4.1 and Satz 4.2 of [TQ there 
exist neighborhoods (5(0) of and (S(w) of uj such that T^iis) is holomorphic in 
5(0) X S{uj). At the same time by Satz 5.3 of [TB], for any oj E Hq the solution 
can be analytically continued to the whole of 6{s{I)). Thus T^{s) is holomorphic in 
d{M.) X Ha. Now restrict back to the real axis and fix the interval TZ. The reasoning 
of the previous proposition can be repeated literally except that now A-jz (a) is not 
zero but equals Ati{uj) = |3a;| sup7j{/^C} < asup7j{/^C}. This results in a similar 
formula as in Proposition 12.41 with perhaps different coefficients, and that proves 
our assertion. □ 

We have an immediate corollary. 

Corollary 2.3. Under the assumptions of Proposition \2.5l if p,q G ^(Ha) then 
for each s £ TZ, rs £ A{Aa)- 
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2.5. Mode decomposition of weak solutions 

The aim of this section wih be to extend the mode decomposition of Solo{T) ob- 
tained previously to entire I?(T)o. Here we assume all the conditions of Proposi- 
tion [53] are satisfied, and we have chosen the system {u^v}a with Ua = TaCa and 
Va = TaCa-, which Span SoloiT). For convenience we will also assume at least the 
part (iv) of the definition of the conventional Fourier transform to hold. 

Unfortunately we do not have a precise analytical description of the Fourier 
transformed test function space even under the assumptions of the conven- 

tional Fourier transform, as it was, for instance, in the Euclidean space by Paley- 
Wiener theorem. In particular we need to know for which modes Ta it holds 

T^{t)~f{a) e V/>) e t e I. (2.19) 

At least we are able to find a sufficient condition under additional assumptions. 

Proposition 2.6. Suppose the bundle T is analytic and Dy,^ has a strictly uniform 
spectrum. For each a e set T'q(O) — p^{uj(a)) and 2^(0) = g'(cj(a)), where 
p\q' e A[Mo]. Then (Eq\2A^ holds. 

Proof. Choose the interval TZ such that it contains both and t. First we note 
that by Corollary EH for a e we have Ta{t) = r\{uj{a)) with r\ e A^q]. Denote 
^^*(A) = rKctt) ) e ^[Ho]. Obviously for any feV{Tt), 

Fi{X^{t))f{a) = [F[{D^,)f]ia), 
where i^/(Dst) is defined by functional calculus. Then by Proposition ! Appendix C.T] 

FiiUt))fia) e P(E). 

Let {[/„} be a covering by local trivializations of Tt, and let be a subor- 
dinate partition of unity. The support of / is covered by Nf (finite) trivializing 
neighborhoods, and we write / = X^n*"/ = J2nfn- It follows / = X^n /« and 
Ta{t)f{a) — J2n '^a{t)fn{oi)- Cousidcr /„ as a section in the trivial bundle 7r~^(?7„). 
As we have seen already (and as we will see even more evidently for homogeneous 
spacetimes in the next chapter) each component supports the Fourier transform 
of one fiber component in some local frame. Thus we can write /„ = where 

G V{Un) and /* is supported in E*. We get 

n i n i 

which completes the proof. □ 

Remark 2.6. An argument involving local trivializations as in the proof of Propo- 
sition [5?n] will show that the multiplication of J-[f{t, .)] by Ia{t) amounts to multi- 
plication of each fiber component by a number, hence Ia{t)F[f{t, .)] e ^'(E) for all 
t e I. 
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Two useful facts about the time dependent Fourier transform can be given by 
the foUowing 

Proposition 2.7. Let f{t,a) e C^^ (l,'D{Y.)j . then 

(i) f{t,x)^F-'[}{t,a)\eV{T) 
(it) J^dtf{t,a) e 2?(S). 

Proof. Let 

f{t,x)=J^~'[f{t,a)] = [ dt,{a)f{t,a)US). 

For each t e I we have f{t, a) £ and therefore f{t, x) E !){%)■ If the compact 

interval A d is such that ^ A, f{t, a) = 0, then obviously Vt ^ A, f{t, x) = 0. 
Because the integration converges in differentiation can be interchanged 

with the integral, thus f{t, x) is smooth in t. The part (i) is proven. 
Now write 

fit,a)= / dfj.ho{x){Ca{x)J{t,x))go, 

and 

/ dtf{t,a)^ dt dfiho{x){Caix),f{t,x))gO = 
Jx Jl JSt 

= / dfj,hoix){Caix), / dtf{t,x))gO=J^[ dtf{t,x)], 

where Fubini's theorem was used with the justification that both integrals run over 
compact supports [T7]- For the part (ii) it remains to show that Jj- dtf{t, x) £ 2?(7t). 
But this is again clear because the integral runs over a compact support. □ 

Next we want to show that the Cauchy problem can be well-posed in the dis- 
tributional sense. We will do it by generalizing fEa l2.ip to distributional solutions. 

Proposition 2.8. For any mo,wi G 'DiTt)' there exists a unique j(mo,wi) — u £ 
V^TYq such that 

u{f) = uoiiU^tEim - miitism, V/ e ViT). 

Proof. By Proposition l2 . 1 1 we know that E is surjective, so we denote the bijective 
part of E to be E^ : I?(T)/ker£' SolniT). For surjectivity of j it suffices to set 

uo{vi) = u{E:^'^ [j{0,vi)]), ui{vo) = -u{E^'^[j{vo,0)]), Vuo,wi e ^{Tt). 

Indeed, 

u{f) = u{E7'[E^[m = u{Er'[j{e,{E^[f]),eAytE^[f]m = 
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= ^i(i^J-lb(^t*(£;J[/]),0)])+^i(i^J-l[XO,^^Vti^J[/]))]) = ^io(^t*(Vti^[/]))-^.l(^*(i?[/])). 

For injectivity of j let uq, ui € ^{Tt)' be given. Define u as in the statement. Then 
obviously u{Df) = because EDf = for any / e V{T), hence u e ViTtYa- Now 
suppose the same formula holds also for different Uq, u'l e I?(7t)' with the same m. 
Then we have 

= K - «[,)(*:(v*i?[/])) - ~ u',)iiUEm, V/ e i?(r). 

Evaluating on / = E^^{j{va,0)) and g = E^^{j{0,vi)) for arbitrary wo,wi G 2?(7t) 
we find uo = Uq and ui = w'j^. □ 

Now we come to the main assertion. Let the modes be chosen such that 
(EqEHl) holds. 

Proposition 2.9. Under the assumptions made, there exist closed topological sub- 
spaces 2?''(S]), 2?^(S) C such that for any ip € 2?(T)o there are unique dis- 
tributions a"^ e 2?"(f;)', b'l' e with 

i^if) = a'^Ka)) + b^Mf)), V/ e 2?(r). 

Proof. Considered as distributions, the functions Ua act as 

uM) = (wa,/)A./ = j^diT„(t)(C„,/(i,.))s, = j^dtT„(t)(C_„,f/(i,.))s, = 

dtT^{t)I^{t)T[rf{t, .)]{-a), V/ e 2?(r). 



(Remember that goo = 1-) The action of Va is similar. By assumption (Ea l2.19|) 
and Remark 12.61 we find 

T^{t)Ut)T[tf{t,.)]{-a)eT>{t). 



Then by Proposition 12 . 71 we get Ua{f) G f (E) (similarly for Va). 

In general, the maps / — > Mq(/) and / Va{f) need not be surjective. Therefore 
we define 

By continuity of the map / — Ua{f) (which is easy to establish), 2?"(S) is a closed 
subspace of 2?(E). Similarly we define I?"(E). 

Recall the mode expansion for arbitrary (p G SoIq{T), 

(j){t,x) = / d^i{a) [aiua{x)+biva{x)\ . 

Thus SoIq{T) can be written as a direct sum of linear subspaces, SoIq{T) = 
Sol^{T)®Soll{T), with 

Soq{r) = {0 e Solo{T): 6"* = 0}, Soll{T) - {0 e 5o/o(r): a"* = 0}, 
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and we will write (p — (p"^ + (p" . Regarding as a distribution in ViJ')' , (j)^ and (j)^ act 
as 

rif) = I dv{a)aiuMl <t>''{f) = I dv{a)bivM)- 

The functions a%,h% can be regarded as distributions a'*' € f)'^{t)', h'*' G P"(E)', 
and we can write 

m = rif) + <l>'{.f) - a'f'iuMj) + b^vM))- (2.20) 

Now let (p e Sol{T) be a solution, which does not necessarily have supp{ip} fl St 
compact. Its Cauchy data are 

{i*{^),i*{Vt'p)) = e S{Tt)®e{Tt). 

Choosing a countable (compactly finite) partition of unity on E we can write 

■DO 

(^0, <^i) = Y^^cPl, 01), {4>l, 01) G V{%) © V{%), 

i=l 

where the sum involves finite items on any compact region f/ G S. Now for each i 
we have 

4,'=^\cl,lcl,{)&Soh{T), 

thus 
Set 

OO OO 

i=l i=l 

and 

OO OO 

vp'' = 5^^^^(0?,■^0r) = $^0^■^ 

i=l i=l 

where the sums converge in compact topology. (This can be seen as follows. The 
intersection of the causal cone of any compact region with a Cauchy surface is a 
compact surface, and therefore only finite summands survive.) But we have 

for some distributions and b^^ . Thus we obtain 

OO OO 

i=l 1=1 

and 

OO OO 

<p''(/) = E'^''''(/) = E^^'(^«(/))- 

i=l i=l 



December 12, 2012 2:19 WSPC/INSTRUCTION FILE Bulk 



30 

This convergence defines distributions 

oo oo 
i=l i=l 

such that 

and thus Sol{T) = Sol'^{T) ® Sol''{T), where 

5'o/"(r) = {(y5 e Sol{T): bf = 0}, S'o/^(r) Sol{T): a'^ = 0}. 

Now let -0 g 2'(T)g be a weak solution, and {xm} a usual moUifier on Sj. Define 
the mollifications ^mV' £ Sol{T) by 

= j(Xm^/'0,XmV'l), 

where V = j('0o,V'i) by Proposition 12.81 Then it is easy to see that '^mip tp in 
V{T)'. That e Sol{T) it follows 

^n.^/- = + (^m^)* G Sol'{T). 

The disjointness Sol''{T)r\SoF { T) = im plies t hat ^ i/-" and ^ 

V'" with some distributions V" G Sol^{T), i/"" e Sol'"{T), such that ?A = -0" + V'"- 
We denote 

0"(/) = hm (9,„V)"(/) = lim a'^-Mf)) = a"/" 

m— >oo m— >-oo 

for some distributions and 6^. Finally we arrive at 

V-C/) = a'"K(/)) + b^{vM)): V/ e i?(r). 

The map 

is a bijection by construction. □ 



2.6. The propagator 

In this section we will find the explicit form of the propagator E in terms of the 
mode decomposition. Of course. Green's functions can be calculated using the tech- 
niques of inverse operators. But our approach will be more concordant to the spirit 
of this work and will at the same time demonstrate the usefulness of the mode 
decomposition in general. 

To use the mode decomposition for weak solutions we assume that at least the 
condition (iv) of the conventional Fourier transform holds, and that the assumptions 
of Proposition 12.91 are satisfied. Choose mode solutions to be such that Tq(0) = 
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T_a(0) and Ta{0) = T_c((0). Then because a — )■ —a preserves both Xa{t) and the 
component E% we have the same mode equations for T„ and T_a, hence everywhere 

T„(t) =T_„(i). 
The function 

detW[T„,f^]{t) = Ut) [f„{t)f„{t) - T„(t)f„(i)] e C°^{I,i- R) 

is the Wronskian of two independent solutions T„ and Ta and is therefore an imag- 
inary constant. For convenience we appoint once and forever to consider only the 
modes normalized by 

T„(i)r„(t) - T„(t)T„(t) = i ■ I^\t). (2.21) 

It can be seen that this condition is consistent with our previous assumptions for 
the modes T„. 

We remark that the Krein space involution f commutes with the connection 
components Fj. Indeed, by definition T — P+ — . where are the projections 
onto the subspaces of positive/negative definiteness of the metric (, )g. Let {e^} be a 
pseudo-orthonormal moving frame, i.e., (ej, ei)g = ±1. The value of each (e^, ei)g is 
preserved under V, and therefore remains in the same eigenspace of F, although 
in our main frame Cj experiences gradient, 

4 n 

j=i k=i 

Hence f commutes with all Fj. We have that 

(w,i')Et = (fu,t;)Et = / dii{a)s{a)Ia{t)u{-a)v{a), Vu e v e ©(7^), 

where s(a) is the Fourier image of the Krein involution F, which due to the remark 
above satisfies s(S') = {+1, —1}, i.e., is constant on each component S*. We have 
used the fact that u(a) = u(— a) which follows from the condition (iv) of the 
conventional Fourier transform. 

Now the propagator is the unique operator E : T>{T) — >■ SoIq{T) which satisfies 

v{f) = {v{t;.),E[m;.)}j:, - {v{t; .), 

Vu e SoloiT), f e V{T), t e I. (2.22) 

As V G Solo{T) we can write 

v{x) = / diJ,{a)a'" {a)ua{x) + / diJ,{a)b'" {a)va{x), (2.23) 

and 

^){a) = a^{a)T^{t) + V{a)f^{t), 
HtTXa) = ay{a)f„{t) + V{a)f^{t). 
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Similarly for E[f] e SoloiT), 



E[f]{x)^ / d^l{a)a^[f]{a)u^ix)+ / d^^{a)b^[f]{a)vo.{x), 



E[f]{t; .){a) = a^[f]ia)T^{t) + [f]{a)f^{t) , 



and 



E[f]{t; .){a) = a^[/](a)T„(t) + 6^[/](a)r„(t), 
with some distribution fields a^[f]{a) and b^[f]{a). Using all this we compute 

{v{t;.),E[fm-)h.-{Ht;-),E[f]{t;-)h. = 



by normalization fEa l2.2l( ) 

- -^ / d^,{a)s{a)Ut) [a"(-a)&^[/](a) - 6''(-a)a^ [/](«)] . (2.24) 



On the other hand, we know that ED ^ 0, thus a^[/](a) and 6^ [/](«) are weak 
solutions of the field equation and can be mode decomposed as 



a''[f]ia) = ai{uf3{f))+al{vp{f)), 
b''[f]{a)=bi{u^{f))+bUv^{f)). 

By (EqI1231) we have 

v{f) = f dfi{a)a''{a)uM) + f dti{a)b'' {a)vM)- 

Inserting (EqU^l), (EqO^B and (EqU^l]) into (EqU^l]) we obtain 

al^bl= 0, al = -bi = i- s{a)6{P - a). 
And our final formula is 



(2.25) 



(2.26) 



-£'[/] (2;) = i / dll{a)s{a) [V-a{f)Ua{x) - U-a{f)Va{x)] 



which is in full accord with the result obtained by ^ for scalar fields on FRW 
spacetimes. 
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3. Aspects of Harmonic Analysis in Homogeneous Spacetimes 
3.1. Spatially homogeneous cosmological models 

The main goal of the current work is to refurbish the mathematical framework of 
quantum field theory on classical cosmological spacetimes, in general, and to ad- 
vance towards a satisfactory rigorous description of cosmological particle creation 
in states of low energy for hyperbolic fields, in particular. The latter would be 
an extension of results obtained in [T] for the Klein-Gordon field on specific FRW 
models to more general situations. Thus although some results were and will be 
obtained under abstract general assumptions, our attention is concentrated at the 
geometrical setup of most common cosmological models. Supported by observations 
of the universe at large scale, cosmology considers mainly spatially homogeneous, 
or in addition also isotropic, spacetimes. A condensed account of cosmological ar- 
guments and their geometrical implications can be found, for instance, in |18| . The 
essence of these geometrical restrictions is mathematically expressed by imposing 
the existence of a sufficiently rich system of symmetries (more precisely, a group of 
spatial isometrics) on the spacetime. Extensive treatments of all possible isometry 
groups and related questions can be found in [18] , [1^ , [2^ . An introduction to the 
generalities of harmonic analysis on vector bundles is given in [21] . In this section 
we will try to deductively introduce our geometrical setup with the help of the 
information in the above mentioned references. 

Foliation by equal time Cauchy hypersurfaces. Recall that we are working 
with a four dimensional globally hyperbolic Lorentzian manifold {M, g) on which a 
global smooth time function and an atlas can be chosen following [8 such that M 
is foliated by three dimensional spacelike equal-time smooth Cauchy hypersurfaces 
and 

ds^ = goodt^ — dcr^ , 

where da^ is the line element on any of those Cauchy surfaces being Riemannian 
submanifolds. 

The structure group. Any vector bundle T can be considered as associated 
to its frame bundle Vj- with structure group GL{n). If we want the fiberwise trans- 
formations to respect the fiber metric, then we have to restrict the principal bundle 
to the orthogonal frame bundle. All fibers Vp with their respective non-degenerate 
pseudo-Riemannian structures Qp are isomorphic, and their generalized orthogonal 
groups 0{qp) (i.e., groups of invertible endomorphisms of Vp preserving Qp) are iso- 
morphic to the generalized Lorentz group 0(±g), where ±g in this context will be 
understood as the signature of q. But the same vector bundle T can be associated 
also to another principal bundle (which we again denote by V-r) with structure 
group H (say, for field theoretical reasons). Then we have a representation r oi H 
on V . If r also respects the metric, then r{H) E 0(±g), so H is homomorphic to 
0(±g). For instance, H = 5*0+ (±g) (tensor bundle) or H — Spin^{±g) (spinor 
bundle) . 
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Isometries. Let us start with reminding some definitions. An isometry of the 
spacetime {M,g) is a diffeomorphism ifj : M ^ AI such that tp*g = g holds on M, 
where tp* is the puUback of -(/'■ If "0' : — > is another isometry, then obviously 
such is also their superposition i/jo?/;'. With the superposition as product, isometries 
thus constitute an abstract group, which we will denote Iso(M). If T — >■ M is a the 
vector bundle over M as defined previously, then an isometry of the vector bundle T 
is a morphism : T — >■ T covering an isometry of the base, tt o v]/ o tt^^ G Iso(M), 
such that 'i>*g — g and \E'*V = V (or more precisely 'i>*D — D when a normal 
hyperbolic field operator D is specified), where ^E** denotes puUback maps, g is 
the pseudo-Riemannian fiber metric, and V is the metric connection. Again via 
superposition, the isometries of the bundle T comprise an abstract group Iso(T). 

The map Iso(T) S^"— )-7ro^'o7r~^€ Iso(M) gives a homomorphism of Iso(T) 
into Iso(Af). The image of this homomorphosm is a subgroup of Iso(A/) and will 
be denoted by Iso^(Af) C Iso(A/), and its kernel is a normal subgroup of Iso(T). 
This kernel Iso(T)/Iso^(Af) consists of isometries of the bundle T covering the 
identity map of M. These are precisely the smooth sections in the principle bundle 
Vr ^ M X H of r, i.e., Iso(r)/Iso'^(A/) = C°°{Vr)- The group muhiplication 
is given by the pointwisc multiplication of sections. 

Homogeneous bundle structure. If the sections in the bundle T are going 
to represent physical fields, than one should have a concrete picture of how they 
transform under the diffeomorphisms of the spacetime M. In case of the tensor 
bundle this picture is automatically encoded in the pullback map. An abstract 
vector bundle does not have such a structure by itself. Thus a physical field theory 
has to specify a homomorphism p : Diff(Af) — > C°°{Vt)- For the tangent bundle 
p{ip) = dip, ij) G Diff(Af). When considering arbitrary diffeomorphisms, then the 
structure group should be GL{n) rather than a smaller H . But if we restrict p to 
p : Iso'^(A^) ^ C°°{Vt), then H can be chosen. For brevity denote G = Iso^(Af). 
We have the injection 

G5g-^gxp{g)eIso{T), 

which gives sense to the left action of G on T by isometries. 

The abstract group of isometries of a pseudo-Riemannian manifold of dimension 
TO is given the compact open topology, in which it becomes a Lie group of dimension 
at most n{n + l)/2 [22]. It can be further shown, that the compact open topology 
in this case is equivalent to the pointwise convergence topology of isometries. Thus 
we automatically obtain a Lie group structure on Iso(Af). Then G C Iso(A^) is a 
topological subgroup defined by 

G^i^e Iso(A/): (V X p{^))* = 0, X pii,))* D = D}. 

If p is a continuous homomorphism, then all the operations in the equations 

X p{iP))* = 0, {ipx p{i}))* D = D 

are continuous, and therefore the subspace G of Iso(Af) defined by this equation 
is a closed topological subspace. But then by Cartan's theorem G is actually a Lie 
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subgroup, as it is a closed topological subgroup ol the Lie group Iso(M). Thus we 
have the structure of a G-homogeneous vector bundle T. 

Spatially homogeneous bundle. The bundle T will be called spatially ho- 
mogeneous if the orbits of Iso^(Af) are three dimensional smooth spacelike hy- 
persurfaces which foliate M . (Maybe it is worth mentioning here that everywhere 
in this work we consider only connected spacetimes M.) By Theorem 8.16 in [2D] 
there exists a parametrization of these orbits by the affine parameter of the family 
of normal geodesies, such that the metric takes the form 

On the other hand, our original foliation by equal time Cauchy surfaces due to 
Theorem 1.1 in |5] also yielded such a metric form. We assume that the time 
function can be chosen such that equal time Cauchy surfaces are the orbits of 
Iso'''(M) (probably this can be shown to be true in general). We note that due to 
the transitive action of G on St for every t, it holds G C Iso'^'(St). We did not 
write G = Iso^*(I]t) because it is possible that for some t ^ t' G 2, Iso^*(I]t) ^ 
Iso"^'' (Sf), i.e., for some time instances the time slice may be more symmetric 
than usual. We will concentrate on G, which is the maximal guaranteed amount of 
symmetry which is present at any time. Thus we see that 7t also has the structure 
of a G-homogeneous vector bundle. 

Consider the principle bundle Vj-t of 7t, which is a subbundle oIVt- The smooth 
left action of G on 7t gives a smooth left action of G on as well. This action 
allows one to construct a global smooth section in , whence it follows that the 
bundle 7* is trivial. Because M ~ Et x I, the whole bundle T is also trivial. Thus 
spatially homogeneous vector bundles over M are necessarily trivial. 

The requirement that the field operator D is G-invariant implies that the func- 
tion m*(x) is in fact a function of time only. 

Homogeneous space structure. Now let StabIso'''(p) C G be the stabilizer 
of G at some fixed point p g M. Then Stablso^(p) is a closed Lie subgroup by Car- 
tan's theorem. That for all p G M, the groups Stablso^(p) are isomorphic, then we 
denote them all by Stablso^(M). In this case the orbits Et of G are diffeomorphic 
to the homogeneous space G/StabIso'^(Af) = S. Denote O = StabIso'^(M)+, 
the identity component. Then F — Stablso^ (M)/0 is a discrete normal subgroup 
of G. If the homogeneous space E is itself a Lie subgroup of G, then it acts on each 
Et simply transitively. 

The four dimensional reality. As already mentioned, the isometry group 
Iso(M) of the n ~ 4 dimensional spacetime M is a Lie group of dimension at most 
n{n + l)/2 = 10. Thus in principle one can construct all real Lie algebras Q of 
dimension up to 10, their corresponding connected simply connected Lie groups G, 
then all discrete normal subgroups F of such G etc., thereby exhausting all possible 
isometry groups of M. This heavy task have been done by Petrov et al [TH] and 
others [50], and all the possibilities are listed in tables. It turned out that only 
the Minkowski space has isometry group of maximal dimension 10, which is the 
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Poincare group. Among all possibilities we are interested in those whose orbits are 
St. Thus the dimension of G is at least three. There are three possibilities of six 
dimensional such isometry groups, which correspond to FRW spacetimes. A number 
of possibilities are available with four dimensional groups, which correspond to the 
LRS spacetimes. And finally there are nine classes of three dimensional real Lie 
groups Bi{N) (called Bianchi groups), which together with their factors Bi{N)/T 
by discrete subgroups F represent the isometry groups of the spatially homogeneous 
spacetimes. It turned out further, that in all these cases besides one (the so called 
Kantowski-Sachs model) the isometry group is the semidirect product G = S x O. 
In this case we will call 7t a semidirect homogeneous vector bundle. In particular, for 
six dimensional FRW groups, four dimensional LRS groups and three dimensional 
Bianchi groups O — 50(3), 5*0(2) and {!}, respectively. The normal subgroups S 
are nothing else than Bi{N)/T. 

In next sections we will work on the semidirect homogeneous vector bundles. 
After establishing the necessary mathematical framework, we will obtain results 
concerning the structure of G-invariant homogeneous bi-distributions. 

3.2. On harmonic analysis in semidirect homogeneous vector 
bundles 

In this section we will collect information on harmonic analysis in G-homogeneous 
vector bundles T GjO where G = E x: O which will be useful later in the work. 
This does not pretend to be self-contained or systematic; quite the contrary, we 
will introduce mainly what we were not able to find in the literature. Otherwise 
references will be provided. 

Semidirect homogeneous vector bundles. Let G = E x O, where O is 
a compact connected type I Lie subgroup, and E a connected normal type I Lie 
subgroup. Moreover, we demand that the modular function of E has a non-trivial 
kernel, so that the representation theories of both E and G are well under control by 
Theorem 7.50 of [23]. We note that this is the case for all Bianchi groups which are 
in fact the only candidates for E in our context. Let E = G/0 have a Riemannian 
structure h which is invariant under the left action of G. Let further T — > E be an 
n-dimensional (real or complex) vector bundle with standard fiber Y and a pseudo- 
Riemannian fiber metric g. Let there be a smooth left action of G on T covering the 
left multiplication of G on the base, such that the fiber metric is invariant under 
that action. Then we will call T a semidirect G-homogeneous vector bundle. If we 
choose an orthonormal frame {Xi} of Tj^E (or {Yi} of T|i), and drag it throughout 
S using the transitive left action of G, we will obtain G-invariant global smooth 
frame {Xi} in T*E (similarly, {Y,} in T). Thus both T*T, and T are trivial bundles. 
Associated to the Riemannian structure h there is a Laplace operator A acting on 
sections / S G°°{T). 

The regular and quasi-regular representations for the line bundle. 
Suppose T from above is a line bundle, n = 1. The left regular representation Lg 
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of G on C(G) acts as 

Because the Rieniannian structure is G-invariant, the metric measure dx is a left 
Haar measure on G, and hence Lg is a unitary representation on L'^{G). 

Now any point x ^ G can be uniquely written as x — x-^xq, where xs T, 
and xo G O. Let rfxs be the metric driven left G-invariant measure on S, and 
dxo the Lebesgue measure on O normalized to \0\ = 1. Then dx = dx^dxo gives 
a left Haar measure on G. Functions / on S are identified with their right O- 
invariant extensions to G, i.e., f{xo) = f{x) = f{xO), for any x E G, o E O. Thus 
G(E) E C{G) (similarly E L?{G),etc.) and we may consider the restriction 

Ug of the left regular representation Lg on G(G) to G(S). Its action will be given 

by 

Ugfix^O) = f{g-'x^O), yx^ e S, .g e G. 

The representation Ug of G is the left quasi-regular representation, and it is nothing 
else but the induced representation IndQl. Note that for O — {1} we simply have 
G = E and = Ug. 

Neither Lg nor Ug need to be irreducible. The central decomposition of Lg is 

Lg ^ dv{'K)Lg{TT), 

JG 

where v{t^) is the Plancherel measure and Lglir) = tt (8) 1 is the primary repre- 
sentation composed of mult{7:, Lg) = dimTf E [l,oo] copies of tt [23]. The central 
decomposition of Ug will be 

Ug = / d^iiTT)UgiTT), 

where Gs C G, dfi is the spectral measure of Ug and for /i-almost all tt, Ug{'K) 
is a multiple of tt (multiplicities mult{T:, Ug) and the measure d/x(7r) need to be 
determined). The corresponding Hilbert space decompositions are 



L\G) = / dv{Ti)H^®Hf, 

and 



G 



L'^iT.) = r d^ji{'K)H{-K), 

JGj: 



where Hin) ^ ® £muit{n,u,) ^ "H^ «) H^. Here C'^^niTrM,) symbolizes some 
Hilbert space of dimension mult^ir, Ug) which is finite or infinite. 

In the following we will deal with Ug keeping in mind that in case G = E 
everything reduces to Lg. 
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The operator H-^. Consider for any n ^ G the bounded operator 

IItt = / don{o). 
Jo 

Then XItt is self adjoint, 

n;= /" do7r(o)* = [ don{o-^) = n^. 
Jo Jo 

Moreover, because O is uniniodular, we have 

7r(o)n^ = 7r(o) / do'nio')^ / d{oo')TT{oo') ^ ^ U^7r{o), e O, 
Jo Jo 

and hence 11^ is a projection, 

Ill= rfo7r(o)n^ = / doU^ = n^. 
Jo Jo 

is a projection onto an invariant subspace of tt\o- Recall the operator D,^ of |23] 
which satisfied D-^Tr{x) = A2 (a;)7r(x)£'7r, for all a; e G. In particular, we find that 
Z?7r7r(o) = 'n:(p)D-n for all o G O, and consequently, D-^Ht^ = H^Dj^- 

The Fourier transform in G/O. The Fourier transform in E = G/O associ- 
ated to Ug is naturally the restriction of that on G associated to Lg ; for /^-almost 
all TT e Gs 

fin) = n{f)D^ e nin). 
For any / e Go(S) and /i-almost all vr e Gs we have 

7r(/) = / dxs / dxo.f{x^O)'K{xT.)T:{xo) = / da;E/(2^sO)7r(xs)n7r. (3.1) 
is Jo 

As usual we have n{Ugf) = n{Lgf) = 'n-{g)n{f) for g e G, f e Go(S). The 
convolution / * /i has the property that if / G Go(G) and h E Go(E) then f * h G 
Go(S). Moreover, it satisfies 7r(/ * h) = 7r(/)7r(/i). 

The case of arbitrary T- Let now dimF = n > 1. The left quasi-regular 
representation of G on C°°{T) acts as 

Recall the G-invariant orthonormal frame {i^ijf^i in T and write any / e C°°{T) 
as f = J2 /'^»- Using that UJY, = Y, we find 

n 

ujf{x)^Y.^grxY,, 

i=l 

where Ug is the left quasi- regular representation of G on G°°{T,). Thus Uj = 
(BnUg, and the harmonic analysis of Uj is the same as that of Ug except that each 
primary subrepresentation of Uj' is the n-fold copy of the corresponding primary 
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subrepresentation of Ug. Making the identification C^{T) 3 f ^ {/'} £ ®„C;f (S) 
we find the Fourier transform of / G (T) to be 

or to say in words, a matrix with n times more columns than that of a scalar 
function. The inverse Fourier transform will be 

n ^ 

fix) L M^)Tr [D^U^TT*{x)f{n)\ x Y,{x). 

i=l "'Gs 

4. On the Fourier transform of distributions 

Here we will collect miscellaneous facts about distributions and their Fourier trans- 
form, which we did not meet in the literature. We continue working with the semidi- 
rect homogeneous vector bundle T with notations established earlier. 

Let P'^(Ge) be the image of D{T) = C^{T) under the harmonic analytical 
Fourier transform /(xs) — >■ /(tt). As we have already seen, /(tt) = ©n/*(7r), hence 
2)'7"(Gs) = ®nI?(S), where ViT) is the image under the Fourier transform of 
C^(S). V'T{Gy) inherits the topology of D{T) via the Fourier transform, and one 
can consider the Fourier transform of distributions D{T)' 3 u ^ u E 'D'~{GT,y 
given by u{f) = u{f). 

The Fourier transform has the remarkable property that it interchanges the local 
and global behaviors. Namely, the local irregularities of a function / are reflected in 
the decay properties of /(tt) at large tt, and conversely, the behavior at infinity of / 
determines the local regularity of /(tt). The precise description of these phenomena 
requires a thorough functional analytical investigation, which we, unfortunately, 
have no possibility to perform here. 

It is widely known that any distribution restricted to a compact region is of finite 
order. In [24] the general structure of distributions of finite order has been found 
for _D(R"). Following a similar pattern we present here a partial generalization of 
that result. By Proposition [Appendix A.Tj let us choose the topology {Xi,2,P) for 
convenience. 

Proposition 4.1. Lei Tr- be an n- dimensional (complex) pseudo-Riemannian vec- 
tor bundle over a connected parallelizable (pseudo-)Riemannian manifold K, and 
let V &e a fiber metric connection. Every u G D{TKy of finite order has a repre- 
sentation 

<f) = Pc^AX^)f)2, V/ e D{T), 

q<k 

where Fa^q E L^{Tk) and the smallest possible such k is the order of u. 
Proof. By our choice 
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Let k be the order of u, i.e., u is continuous in ||.||fc-norin. Define tlie following linear 
injective map 

V : DiTK) ^ ^ = ^ L^Tk) 

q<k 

by 

V(/)=0Pa,,(X,)/. 

q<k 

Then obviously ||V(/)||<i. = ||/||fc. If we denote by ^' = V {D{Tk)) C then 
uoV^^ is a continuous functional on ^ with the norm ||.||$, and thus by Hahn- 
Banach theorem can be extended to a continuous functional F € But $ is a 
Hilbert space, thus $' = $ and F G and for any G 

g<fe 

This yields our desired formula 

q<k 

If such a formula held for a smaller k, then obviously the order of u would be 
smaller. □ 

Several variations of this proposition may be established by choosing different 
norms. Note that the order of a distribution, if finite, depends on the choice of the 
family of norms defining the topology. 

Remark 4.1. As already mentioned, any distribution is locally of finite order, 
hence the proposition applies to the restriction uk G C^{T\k)' of any u G D{T)' 
to arbitrary compact connected region C S. 

We come back to our homogeneous bundle T and proceed to the Fourier de- 
scription of distributions u G D{T)' of finite order, which again can be applied for 
restrictions to compact regions. 

Proposition 4.2. Any distribution u G D[T)' of finite order is given by 

u{f) = / dis{^)Tr [u(7r)*/>)l , 

where u{-k) : C™"'*^'^'^'')*" -> "Htt is a fi-locally integrable field of Hilbert- Schmidt 
operators. (Note that the trace operator includes also the summation by fiber indices 
i = 1, which now enumerate blocks of columns.) 

Proof. Let k be the order of u. Choose {Xi} to be the generators of left translations 
on C°°{T) and let by Proposition 14. II write u as 

q<k 
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Consider the Fourier transform 

xTm = I d.^ (lim .(x.)n.i?. 

where is the corresponding element of the Lie algebra of E. The integral runs 
over a compact region, and is therefore uniformly absolutely convergent with the 
Hilbert-Schmidt norm, thus we can interchange the limit with the integral, 

-'^i/(7r) = linij / rfa;E(C/oxp(-ie.) - l)/(a;s)7r(xs)n^-D7r = 

^j.^^(exp(-iMl^/(,). 

t^o t ' 

On the right hand we see nothing else but the generator of the derived representation 

of TT, 

y 7r(exp(~te,)) - 1 

hm = ^OjTT, 

t^Q t 



whence we find 



As a result we have 



and thereby 



XJ{7t) = -5,7r/(7r). 



q<k 



d^{'K)Tr 'u(7r)*/(7r) 



where 



i(7r) =^[P„,,(-a,;^)]*F„,,(7r). 



This completes the proof. 



□ 



Such a result should not be surprising. If the measurable functions Fa,q were q 
times differentiable within the space of locally integrable functions, then we could 
hypothetically use integration by parts to make all the terms in the formula of 
Proposition 14. II of order 0, which would correspond to a regular distribution. The 
failure of the derivatives of F^^q to remain locally integrable is reflected in the fact, 
that multiplication of Fa^q{'K) by diii* makes it not square integrable any more, but 
possibly only locally integrable. This reflects the local-to-global interchange made 
by the Fourier transform: higher frequencies feel local irregularities. 

The image 2)^(Gs) of compactly supported smooth sections under the Fourier 
transform is of considerable interest. In harmonic analysis it is described by various 
Paley- Wiener type theorems. Although there are refined Paley- Wiener theorems 
for adapted Fourier transforms for certain classes of semisimple or solvable groups. 
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there seems to be no such theory for the general abstract setup. Next we present 
a partial answer to the problem, namely, a criterion for smoothness for sufficiently 
decaying functions, which gives hints about how the general solution might look 



Proposition 4.3. For a function f £ i^(T) the following two conditions are equiv- 
alent: 

(i) for any polynomial P{Xi) of generators {Xi\ with constant coefficients, 



We can go a step further and establish a weaker necessary condition for a distri- 
bution to be given by a smooth integral kernel. For this purpose we want to remind 
a few definitions on a more abstract level. 

Let be a test function space. We have T>{S) C L^{S) and therefore 

L°°{Sy C ^{Sy . Let {rji} be a finite system of linear maps r/i : S ^ S. A distribu- 
tion u S 'D{Sy is of rapid decay in {rji} if for any polynomial P{r]i) of variables {r]i} 
it holds u{P{r]i).) e L°°{Sy. We will symbolically write this as w = o{{rii}~°°). If u 
is given by a locally integrable kernel, and {rji} are coordinate operators, then this 
definition obviously reduces to the usual criterion for functions of rapid decay. 

Proposition 4.4. For a distribution u € 'D{Ty from V'^ {G^y B u = o{{din}~°°) 
it follows that u has a smooth integral kernel. 

Proof. That u is smooth means that all derivatives of all fiber components 
are continuous. In other words, for any polynomial P in the generators {^j}, the 
distributions P{Xi)u^ can be evaluated pointwise. A precise statement can be given 
as follows, u is smooth if and only if for any polynomial P{Xi), point tn E S and 
sequence of test fimctions fg S{x — m) in C5°(E)', the following limit exists for 
all j — 1, n and is finite, limg_).oo {P{—Xi)fq). The Fourier transform of the 
distribution 5m = 5{x — m) can be easily read from the Fourier inversion formula. 



like. 




Proof. As we have seen in the proof of the previous proposition. 





That means fq 5{x — m) is equivalent to fq — )• 7r(m)n^Z)^ in the weak sense. 
Hence 



P{-Xi)fq ^ P{diTT)TT{m)U^D^ 
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in the weak topology. It follows 

lim u'{Pi-Xi)fq) = {P{d,n)7r{m)U^D,) (4.1) 

f oo 

whenever one of the sides converges. 

Now suppose u = o({5i7r}-°°). Then for any / e £°°(Gs) (i.e., ||/(7r)|| e 
L°°{Gs) in the usual sense) we have 

M^' (P(a,^)/^(7r)) <oo, j-l,...,n. 
In particular, 71(171)11^^0,^ G L°°{Gj:), whence fEa l4.ip follows. □ 

We are incline to think that this necessary condition is not far from the desirable 
equivalent condition. This is, however, an open problem in harmonic analysis, and 
we only hope to be able to give a satisfactory answer in the future at least in the 
context we are interested in. 

5. The adapted Fourier transform 

We start by noting that because the function m*{t) is a function of time only, 
the eigenfunctions of are nothing else but the eigenfunctions of the Laplace 
operator Aj. In the first chapter we introduced the eigenfunction decomposition 
associated to any self adjoint operator as the Laplace operator A, 

/^/(a) = Ca(/), 

where Cq-s are the generalized eigenfunctions of A. Putting additional structure 
related with particular geometries one arrives at various Fourier transforms, which 
are very practical in many respects. On the other hand, the abstract harmonic 
analytical Fourier transform is a powerful tool for analyzing general problems and 
properties, but its machinery is functional analytically complicated for use. These 
two theories are, however, related, although the exact relations have not been suffi- 
ciently explored in the literature so far except for compact groups. In the compact 
case the eigenfunctions of A are the matrix elements of the irreducible represen- 
tations for some choice of the basis, and the two techniques can be unified. Each 
choice of the basis results in a Fourier transform which is adapted to it, hence such 
transforms are sometimes called adapted Fourier transforms. In the non-compact 
case functional analytical complications arise, though intuitively the situation re- 
mains similar. In this section we will try to construct adapted Fourier transforms 
at least on our semidirect homogeneous bundle T. 

The Laplace operator A is invariant under G and hence commutes with . 
This means on each primary component it acts as a multiplication from the right 
by a possibly unbounded self-adjoint operator A(7r), 

A/(7r) = />)A(^). 
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For any / G L^{T) we have that A/ is a distribution of order at most two. By 
Proposition 14.21 it means that the muhipUcation of any Hilbert-Schmidt operator 
/(tt) by A(7r) from the right leaves it again Hilbert-Schmidt. Let (t{tt) C M be the 
spectrum of the self- adjoint operator A(7r) as acting from the right (this spectrum 
is non-positive, because A is an elliptic operator). For each A G cr(7r) let ^7r,A,r,s be 
the generalized eigenfunctions of A(7r), i.e., distributions satisfying ^7r,A,r,s^(''i') = 
A^7r,A,r.s which are linearly independent and complete in H(7r) for r G R-^ C M and 
s d S!^ <Z M. (they can be constructed from delta functions using the spectral 
theorem). Now consider the following distributions in the Fourier space. 

Their preimages are distributions C-n-.A.r.s G D{T)' which are generalized eigenfunc- 
tions of A, and by elliptic regularity theorem, are smooth sections in T. Thus we 
have found, that the adapted Fourier transform /(tt. A, r, s) is nothing else but the 
coefficients of /(tt) as expended in the system ^7r,A.r,s- It is worth noting that r 
parameterizes "Htt, and A, s parameterize C™"'U'r:C^s) ,|, Actually, S^^ consists of 
n copies of some set S^^a- 

The choice of the system ^7r,A,r,s is rather arbitrary and leaves room for adapta- 
tions. The first adaptation we wish to make is the following. For any Cir.A.r.s we want 
C,T,\,r.s — C-k' ,\' .r' ,s' for somc other parameters. Obviously A = A', and it is easy to 
see from the Fourier inversion formula, that this amounts to requiring that ^7r,A,r,s 
enters the system ^7f,A,r',s' for the representation tt with some other parameters r', 
s' . The representation tt may lie in the same equivalence class [tt] or not. 

Lie groups are analytic manifolds, and all the group and algebra structure is 
given by analytic functions in any analytic atlas. In particular, the eigenfunction 
problem /S.C,T^^\_r,s = ^C-n.\,r.s is an analytic elliptic equation, and the solutions 
C7r,A,r,s(a^) are therefore analytic functions in x. If S is compact, then Gs is discrete, 
and each (T(7r) is also discrete. Representations are finite dimensional, hence r and 
s run over finite sets. The set S = {a = (tt, A,r, s)} can be considered a discrete 
manifold symbolically divided into n components as corresponding to each copy of 
S'tt^a- The space 'D^{G^) corresponds now to the space I?(S) of functions on E, 
which are of rapid decay in A, and also decay sufficiently fast in tt by Proposition 
1431 

If E is non-compact, suppose there exists a subset K C Gs such that 
/x(Gs \ K) — and K can be cast into an analytic manifold. Then we can restrict 
our Fourier transform from Gs to K without violation of the Plancherel equality. 
Suppose further that the set — {a — (tt. A, r, s)} can be made an analytic man- 
ifold consisting of n disjoint components as in the compact case. Each component 
itself may have several connected components if 1 < •mult{'K, Ug) < oo, in which 
case s will run over a discrete set. Then we can choose Cir, a, r.s to be analytic in 
all its parameters (if s is discrete, analyticity in s is void), so that V'^ {G^,) will 
correspond to the space £'(E) of some analytic functions on E which have at least 
above mentioned decay properties in A and tt, but also are in r, and in s if the 
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latter is continuous. 

Finally let us define a symbolic involution a — > —a on E satisfying = 
^Q. Clearly this involution will preserve A. Now if the necessary assumptions are 
satisfied, we arrive at a conventional Fourier transform. In the next publication we 
will see that in the majority of situations in cosmology these assumptions are valid, 
and that will enable us to exploit the machinery of mode decomposition to our 
cosmological models. 

5.1. Invariant bi- distributions 

In this section we will try to analyze the structure of bi-distributions w £ 
{T>[T) (8" ViT))' which are invariant under the left quasi-regular action Uj of G on 

w{Ujf, Ujh) = w{f, h), V/, h € V{T), 

and compare with results obtained earlier in the literature. 

Decomposing each / = Y^P^i, p € we find for u e V{T)' and w G 

(2?(r) ® V{T))' 

71 n 

e Co"(S)', e (Co"(S) <E> Co°°(I]))' , 

so that the problem reduces to that for scalar distributions. 

The following proposition establishes the general form of the G-invariant (or 
homogeneous) bi-distributions. Our approach is greatly inspired by ^ where this 
analysis is performed for M". 

Proposition 5.1. Every w E (C^(S) (g) C(f(E))' satisfying w{Ugf,Ugh) = 
h), V/, h e C^{T.), geG, has the form 

w{f, h) = Uy, (/* * h) 

for some G C5"(I])'. And conversely, any G C^CE)' gives rise to such an 
invariant bi- distribution w. 

Proof. Recall that for scalar functions Ugf{xY.O) = f{g~^x^O). By the nuclear 
theorem w can be uniquely extended to w e C^^E x S)' via embedding 

C^iE) ® Co°°(E) 3 /(xs) ^ h{y^) ^ f{x^)h{y^) e C^{G x G). 

That 

w{f{g-\^0),h{g-^y^O)) - w{f, h) 
by continuity implies that 

w[c^{g-^x^O,g-^y^O)) = w{c^{x^,y^)), V0 G C^[E x E). 
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Define the linear automorphism 

by 

Mxi:,yi:)= / dxo(p{xY„XT.xoyT.O) = / dxo(f>{xj:,x^Oxoyi:0), 
Jo Jo 

and the puUback of w under this automorphism by i), ^(V'^) = w{4>). If <pg{x-E,yj:) = 
^{g~^xj:0,g~^yj:0) then 

V'0s(a;E,yE) = / dxo(l){g~^xj:0,g~^xj:Oxoyi:0) =-ilj^{g~^x-EO,yj:). 
Jo 

Now 
thus 

Consider the restrietion f of w to C^(S) ® C(y^(S). The last equation implies 
v{f{g-^xj:0),hiy^)) = v{f,h), yf,hG Cg^{^). If we fix h, then v{.,h) e Co~(S)' 
is a distribution which is invariant under all translations, and is thus given by a 
constant kernel, v{f,h) = Uw{h) Jj,dxj:f{x^), for some Uyj : C^(S) — >■ C. On 
the other hand, if we fix /, then continuity in h implies G C^(S)'. Because 
the integral J-^ dxsf{x^) runs over a compact region, it can be transferred into 
Uw, i.e., v{f, h) = Uw {Jy, dxj:f{xj:)h{yj:)). This in turn implies by continuity, that 
v{ip{xj:,ys)) = Uw {J-Ydx-Eip{x-E,yj:)). Finally we arrive at 

w{f,h) = w(/(xs)/i(2ys)) = v{fixi:) / dxoHx^Xoy^O)) = 

Jo 

dxY.f{xT.) / dxoh{xY,xoyT,0) ) = Uw{r * h). 
T. Jo J 

The converse statement is obvious. □ 
For a distribution w G {p{T) X>(T))' this will mean 

n 

Remark 5.1. Note that every G-invariant bi-distribution w G (G^(S) ® G^(S))' 
is in particular S-invariant. Let f * h {f * h) and /* (/*) denote the convolution 
and the involution with respect to G (S), respectively. Then 

w{f,h) = Uni i^j dxT.f{xj:) J dxoh{x^xoy^O)j = 
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l^dxoL^-J*^h{ysO)^ =u'M**h) 

for some other u'^ £ (C^(S) ® C(f (S)) as expected. 

Let 2?(Gs) be the image of C^(E) under the harmonic analytical Fourier trans- 
form /(xe) — /(tt). As an obvious corollary we arrive at the form of an invariant 
bi-distribution in the Fourier space. 

Corollary 5.1. A G-invariant hi- distribution w S {C^ (T,) <Si C^^ (T,))' in the 
Fourier space is given by 

w{f,h) = uMfik^)) = uMfrhiTr)) 

for some iiw G V^G^)' . 

An immediate consequence of Proposition 14. 21 is the following 



Corollary 5.2. Under the assumptions of Proposition \4-^ a G-invariant bi- 
distribution wk G {C^{K) ® C^{K))' with K CT, compact is given by 

WK{f,h)=^ / dfi{T:)Tr (7r)*7r(/)*ft,(7r) 

Proof. It suffices to note that 

supp{/ *h} C 0(supp{/})"^supp{/i}0, 
and to apply Proposition 14.21 □ 

Finally we establish a generalization of the results by 5 for FRW spacetimes. 

Proposition 5.2. Suppose that the group G is such that all multiplicities 
m,ult{Tr, Ug) are finite. Then any G-invariant bi-distribution w G {D{T) ® D{T))' 
has the form 

w{f,h)^ I d^l{n)Tr\{f{^)u{^)yh{^l) 

where u(Tr) is a fi-locally measurable field of [mult(Tr,Ug) ■ n] x [mult(TT,Ug) ■ n] 
complex matrices. 

Proof. Let start with the case w G (C^(S) (8) Cg"(S))'. The condition that the 
modular function of E has a nontrivial kernel ensures that the formula (7.49) of 
(Folland) is valid, so that for //-almost all tt the operator Dj^ is invertible (injective). 
Therefore we can write 7r(/) = f{'K)D:^^, so that 7r(/)*/i(7r) = D:^^ f (tt)* h{Tr) where 
/(7r)*/i(7r) is a mult(n,Ug) x mult{TT,Ug) complex matrix. Now for any compact 
C E by Corollary [Owe find that the restriction of w; to C§°{K) (g) C^{K) 
is given by 

WKif,h)= f d^li7^)Tr\u'J,i7^yD;;'fi^^)*h{^^)\ ^ f d/i^Tr [uif(7r)7>)*/i(7r) 



Ge 



Gy 
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where uk{t^) is a mult{Tr,Ug) x mult{TT,Ug) complex matrix. Choosing a larger 
compact K d K' d we will arrive at another matrix uk'{t^)- But when restricted 
to K, wk' must coincide with wk, hence uk'(t^) = 'U_fc(7r). Thus the matrix uk(t^) 
is the same for any K, and the formula holds for the entire w. 
Now for w e (D(T) (g) D(T))' we have 



which completes the proof. 



dii{ii)Tr (/(7r)u(7r))*/i(7r) 



□ 



In the case of FRW spacetimes all the assumptions of the last proposition are 
satisfied. In particular all mult{TT, Ug) = 1 and for the scalar case we find that 
any G-invariant bi-distribution is given by a locally measurable scalar field u{t:). 
The additional constraint of [5] of being polynomially bounded is an artifact of the 
imposed additional continuity requirement. 
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Appendix A. Space structures. Distributions 

Let us start with introducing symmetric metric products 

(/, h}M - / dfig{x){f{x),h{x))g, f e £{T), h G V{T), 
{f,h)s,= f d^l^{x){f{x),h{x))g, f e£{rt),hev{rt). 

JT. 

The pseudo-Riemannian metric induces a Krein space structure on V, the 
typical fiber of T. Whence there is a Krein involution F, such that (u, w)g = (m, rw)g, 
u,v & V , \s & positive definite hermitian inner product. This gives rise to positive 
definite hermitian inner products 

{f,h)M^ I dfig{x){f{x),h{x))„fe£{r),heV{r), 

JM 

{f,hh,^ [ dfi^{x){f{x),h{x))„ f he ViTt). 
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The completion of spaces I'(T) and with respect to these products becomes 

the Hilbert spaces L^{T) and L'^{Tt), respectively. The tangent space TpM at a 
point p € M with the Lorentzian metric g is another example of a Krein space. In 
the same spirit one defines the positive definite inner product (, )g on TpM. The 
metric h on TpE is Riemannian, so the construction of (, )h is straightforward. Note 
that and g together give pseudo- Riemannian metrics on all product bundles 
T*M <Si ... «) T*M ® r (respectively, (, )g and h on T*E (g) ... <Si T*E (g) 77). All the 
resulting standard fibers are again Krein spaces, and can be given inner products 
(, )g in the same fashion. These in their turn produce products (,)a/ and (,)st on 
the respective sections. 

The perfect countably Banach topology of the test function spaces 'D{T) and 
D{Tt) can be given as usual (e.g., 7 ). However, as we are going to perform a spec- 
tral analysis, we will need nuclear countably Hilbert space structure, to which we 
proceed Let O C M be a compact region. Let 

Vo{T) = {fe V{T): supp{f} c O) 

and define the family of positive definite inner products {,)o.p on DoiT) by 

(/, h)o,p = 5]((V)V, (V)«/i)m, V/, h e Vo{T), p, g e N, 

which induces a family of norms ||.||c),p. One can show that this family of norms is 
growing and consistent, and gives the same topology as the usual one. Let us give 
T^O (T) a countably Hilbert space structure in the following sense, 

N 

It can be shown, that thus constructed countably Hilbert space VoiT) is nuclear. 
Let now 

d C ... C On d ...dM 

be an infinite family of growing compact regions. Then give 'D{T) the inductive 
limit topology 

V{T) = lim I?o„(r). 

Here we are done. Distributions 2?(T)' and operations on them can be defined as 
usual. The same construction can be done for I'(7t) with minor modifications. 

At the end let us consider the choice of the topology in detail. In the literature 
one usually chooses the family of norms ||.||p (or sometimes a family of seminorms 
|(.)|p; from these seminorms one can make norms by |j.||p — X]q<p KOI? o^' II -lip — 
sup^^p |(.)|g etc.) rather arbitrarily in accordance with the setup of the problem, 
and it is tacitly assumed but not everywhere proven, that all such choices give 
equivalent topologies. Let us for consistency present here a proof of this fact. The 
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zest of the proof (the usage of the Sobolev embedding theorem) was suggested by 
G. FoUand. 

Proposition Appendix A.l. Let T M he an n dimensional pseudo- 
Riemannian vector bundle over the d-dimensional parallelizable pseudo-Riemannian 

m,anifoM M with positive metric product (, )g constructed as above, so that we have 
well defined L™ norms |||-|||rn for 1 < m < oo on 2?(T). Let V be a connection on 
T. Let 

(i) X1...X4 be a system of first order smooth differential operators on C°°{T) 
which span the tangent space T*M everywhere 

(ii) the seminorms he given by \{f)\a.q = |||^Q,g(-'^i)/|||m; where Pa,q{Xi) are vari- 
ous m,onom,ia,ls of order q in {Xi}, f € T^{T) 

(Hi) the family of norms he given as \\f\\p — |||{|(/)|Q,q}q<p|||i'=; or by a superposition 
of different 1 < fc < 00. 

Then the topology of'D{T) defined by this family of norms is independent of the 
decisions (i) to (Hi). 

Proof. For convenience denote by (Xj,m, *) the triple of choices at points (i),(ii) 
and (iii). Then (Xj,m, *) (X-,m',*') will mean that this two topologies are 
equivalent. 

As the topology of X'(T) is the inductive limit of various ^{Tk) with Tk = 
■JT~^{K), K C M compact, it suffices to prove the assertion for an arbitrary T>{Tk)- 
The topologies given by two families of norms {|| -Hp} and {|| -Hp} are equivalent if and 
only if these two systems of norms are themselves equivalent, i.e., Vp, 3q{p), r{p) > 0, 
< C7p,C; G R such that < and ||.||; < C'^W.^py Let us start with 

the point (iii). Suppose the choices (i) and (ii) are fixed, i.e., consider (Xi,m, *) 
and (Xi,m, *'). Then all possible choices in (iii) give equivalent systems of norms 
because of the elementary inequalities 

|||{|(/)|a,J/|||;oo < ... < |||{|(/)U,J/||b. < ... < |||{|(/)|a,,}/||bi < ^/|||{|(/)|a,,}/|||;oo, 

where Ni is the number of terms in the index set /. These inequalities can be applied 
consecutively to estimate any composite norm by, say, |||.|||;oo. An example of a 
composite norm is ||/||p = supg<p IHV/llloo- We found that (Ai,m,*) ~ (Aj,m,*'). 

Now let 1 < m < 00 at (ii) and A; = 00 at (iii) be chosen, and choose two 
systems of operators {Xi\ and {Yi} at point (i) to construct the families of norms 
{II.IIj,} and respectively. This corresponds to {Xi,m,l°°) and {Yi,m,l°°). 

Because {Aj} spans T*M, there are functions Cy (x) e C°°(M) and smooth fields 
of homomorphisms Ti € C°° {Hom{T,T)) with Yi{x) = Cij{x)Xj{x)+ri. Using 
this for any monomial Pa,q{Yi) we get 

PaA^)f = J2<,MAXi)f' 

/3 
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where ^(x) G C°°(M) and qi^i) are monomials of order less or equal q. The 
number of summands is less than, say, (4c?)''. It follows by Minkowsky inequality 

and then by Holder inequality 

\\\ci,qix)Qi^g{Xi)f\\\rn < Ca,q^ \\\Qi,qiXi) f\\\m = C'a.q^ \if}a{a,q,l3),q{a,q,l3), 

where < Ca.q = sup^ ll|cQ^q|||oo- In other words, the seminorms of order q of the 
second system can be estimated by linear combinations of seminorms of the first 
system of the same or lower order. Then 

II /lip = SUp|(/)|'„g < Cp sup J2\if')ic,{a,q,P)Mo',q,P) < 

q<p q<p jj 

< Cp(4d)''sup|(/)U(„,,,^),,(„,,,^) < C,(4d)''sup|(/)U,, = Cp(4df ll/llp, 

q<p q<p 

where < Cp = supg<p Ca,q- For the other direction of the estimate we simply need 
to switch {Xi} and {Yi}. Thus these two topologies are equivalent, {Xi,m,l°°) ^ 

(y„m,;-). 

Finally let Xi — Vi (components with respect to a global orthonromal frame 
in T*M) be chosen at (i), and ||.|| — |||{|(-)lQ,g}9<pll|p at (iii). We construct two 
families of norms by choosing 1 < m < oo and to' = oo at (ii) for ||.|jp and ||.||p, 
respectively. This can be symbolized as (Vi,TO, and (Vi,oo,/^). Because K is 
compact, by an application of Holder inequality we obtain 

1 1 1 ■ 1 1 I m ^ C-fyi III . 1 1 I oo 

for some < Cm G M, and hence obviously 

ii-iip <anii.|i;,peNo. 

The opposite inequality requires an application of Sobolev embedding theorem for 
compact manifolds '5S|,'55|. Denote the Sobolev norms (which are equivalent to 
those in [2 5) ) 

/^lllvvillL 

V q<p 

Then an application of Sobolev embedding theorem gives 

|||-|||vi/0-°° = IIMIloo < -Dlll-lllvyd,! 

for some < £> G K. By another application of Holder inequality we find 

||Mlln/'i.i < IIMIIvF''.^, 
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and therefore 



\oo<Dj^\\\V^f\\\l 

q<d 



Next 

ll|vvilli = Elll^""^(^^)/llli' 

a 

and finally 



lip = jEIII^«.'^(^^)/III- ^ ^JEE III v'^«.«(^^)/lll2 

V 9<P V <l<Pi<d, 



y q<pj<d y q<P+d 



<DE \\\p.,q{xm?m = m\f\Ud, 

Y q<p+d 

where in the last inequality again Holders inequality was used with some < € M. 
Thus we have shown that choosing any 1 < m < oo is equivalent to choosing m = oo 
at point (ii), i.e., (Vi,m, Z^) ~ (Vi,oo,Z^). 
Write 

{Xi,m,*) ~ (Xi,m,/°°) ~ (Vi,m,Z°°) ~ {Vi,m,l^) ~ (Vi,oo,/2) ^ {Vi,m',l^) ~ 

~ (Vi,m',/°°) ~ {Xi,m',n ~ {Xl,m',*'). 
The proof is complete. □ 



Appendix B. On the time dependent harmonic oscillator 

Here we will concentrate on some properties of the solutions of the smooth complex 
time dependent harmonic oscillator equation 

f(s)+A(s)T(s) = (B.l) 

where A(s) is a smooth complex function on the real line. This equation is under 
attention since a long time, but some results are not that easily available today (at 
least for us). 

We start with an easy remark. Denote by 
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the Wronski matrix of two solutions Q and R. 

Remark Appendix B.l. Let Q,R be two linearly independent solutions of 
(Ea lB.ip . and T an arbitrary solution. Then from the conservation of detM^[(3,T] 
and det T] it is easy to find 

^ ^[^'^]"'(^) ^ w[Q,Rm X (^^^^^^^ = 

- aetW[Q,R] (0) (^_^^^^ ^^^^ j X ^^^^^ ^^^^j X . 

Thus having at hand two such particular solutions Q,R, we have a control over 
arbitrary solutions T in terms of their initial data. 

Our first task is to obtain a control over the magnitude of the solution T on a 
given compact interval TZ in terms of its initial data T(0) and T(0). This is done 
by the so called energy estimate. Define the energy of a solution T by 

W[f]{s)^^\f\^is) + ^mis)\Tis)\\ 

If 3fiA > on n then 2W[f] dominates 3?A|Tp and |rp, and obtaining bounds on 
W[T] we automatically get bounds on |T| and \T\. 

Proposition Appendix B.l. For arbitrary solution T of 

with smooth complex valued A(s) having a positive real part (i.e., 5RA(s) > 0) on a 
compact intervalTZ, the energy function W[T](s) satisfies the estimate 

W[r](0)e ^ V^^ < W[T]{s) < W[T]{0)e° ^^Ssam^'" 

for all s CzTZ. 

Proof. Write T{s) — R{s) + iS{s), A(s) — 6(s) + iS(s), and insert into the equa- 
tion. We will get the following system of real equations, 

R{s) + e(s)i?(s) - E{s)S{s) = 0, 
Sis) + e{s)Sis) + E{s)R{s) = 0. 

We can cast this into a real vector equation 

f(s) + A(s)f(s) = 

by denoting 

fis)^{R{s),S{s)V, 
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and 



where A='= denote the symmetric and antisymmetric parts. The energy function 
equals 

mT]{s) = lf\s) + iT^(.)A(.)f (.) = ^f\s) + if ^(s)A+(s)f (s). 

On the interval TZ we have W[T'](s) > as by the assumption Q{s) > 0. One can 
easily find that 

VV[f](s) = f ^(s)A-(s)f (s) + if ^(a)A+(s)f (.), 
whence it follows 

W[f]{s) < \E{s)\\f{s)\\f{s)\ + |a,lne(s)|W[f](5). 



< 



By definition of W[f] and positivity of 6 we have |f (s)| < ^2>V[f](s) and |f (s) 
^J2yV[f]{s)/Q{s) on TZ. It follows then 

and integrating this we finally arrive at 

- fn rfo-(^'f' +\ds Ine(o-)l) ^ \ f' drrCf^"''^ +\ds In e(cr)\)\ 

precisely as in the statement. □ 

If however A is not guaranteed to be positive;, then on those regions where it 
is negative the magnitude of the solutions is expected to behave exponentially. We 
are able to capture that exponential factor by the following beautiful trick. 

Proposition Appendix B.2. For any < k e M, any solution of the equation 

f (s) + A{s)T{s) = 

can be represented as T{s) = th(Ks)) ch(/ts), where t{z) is a solution of the 
equation 

f{z) + Q{z)t{z) = 

with 
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Proof. The proof is elementary once we already know the clue: the substitution 
of variables nz = th(Ks). The substitution r(s) = r(s) ch(Ks) into the original 
equation gives 

f(s) + 2Kth(Ks)7-(s) + [k^ + K{s))t{s) = 0, 
then the substitution s — > z yields the final formulas. □ 

Let us say a couple of words about this. If SRA has a minimal negative value 
— c on some domain, then it suffices to set k — ^/c to reduce the problem to an 
oscillatory equation for p. The upper bound of the rate of exponential expansion is 
precisely given by the square root of the minimal negative value of SRA. 

Finally we combine these two statements to find an explicit uniform bound on 
an arbitrary solution T. Let the compact interval TZ containing be fixed, and set 



A-ji = sup |5A|, CTz = inf 5RA, k = yj\ + | niin{0, ck}|, = sup \ds In (k^ + 5RA) I , 

Du = sup(k^ + 3?A), e-Ti = inf(K^ + 5RA) = 1 + max{0, ck}, 
n K 

Ln ^ {2ATz + ch^{n\n\)Bn +2Ksh(2K|7^|)) 
(we suppressed the index TZ oi n for convenience). 
Corollary Appendix B.l. For an arbitrary solution T it holds 



\T{s)\ < |r(0)K/— e^'^ ch(«|7^|) + |T(0)|-Le^^ ch(/^|7^|) 



for all s (zTZ. 

Proof. Consider the linearly independent solutions Q and R given by initial data 

Q(0) = 1, g(0) = 0, i?(0) = 0, R{0) ^ 1. 

Using Proposition [Appendix B.^ represent them as Q{s) = £,{z{s)) ch{Ks) and 
R{s) = p{z{s)) ch(Ks), where ^(z) and p{z) are solutions of the equation 

f{z) + n{z)T{z) =0 



with 



Using z(0) =0 and 



n{z) = 



K + A(i ath(Kz)) 



(l-K^Z^)^ 



^ [r{z{s)) chins)] = + r{z{s)) sh(«s)«, (B.2) 

we find 

e(0) = l,e(0)-0, p(0)-0, p(0) = l. 
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Note that 3firi(s) — k'^ + diA > 1, thus Proposition | Appendix B.Tj is appUcable for 
C and p. We have >V[C](0) = ^^{0) and W[p](0) = i. Now 

dZ as ria 



dz 



ds 



ch^fKs)— In (k^ + ^A(s)) +2Ksh(2Ks). 
ds 



Then it follows 



y'Ml{a) 



< -th{K\n\) {2Atz + ch^{K\n\)BTz + 2Ksh{2K\n\)) 



< 2 {2ATz + ch'^{K\n\)Bn +2Ksh(2K|7^|)) = 22.^. 
By Proposition [Appendix B.I] we have 

W[C]{z) < W[e](0)e2-^«, W[p]{z) < W[p](0)e2-^^, 

which entails 



y/K^ +m{s{z)) 



e^", IpWI <e^ 



For Q and R we get 



|Q(s)| < J— e^-ch(K|7^|), < -^e^-ch(/t|7^|), 



en 



1 

en 



and using fEa lB.2 



\Q{s)\ < ^e-- ( 1 + ^^hmi) , |^(,)| < f 1 + 



Finally let T be an arbitrary solution of the original equation. Applying Remark 
[Appendix B.l for Q,R and T we find 



T{s)^T{0)Q{s)+f{0)R{s), 



and hence 



\T{s)\ < |^(0)u/:^e^'^ch(K|7^|) + |T(0)|-^e^'^ch(/^|7^|), 



en 



1 



as asserted. 



□ 
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Appendix C. A result from functional calculus 

In this section we will obtain a result using the theory of holomorphic functional 
calculus of strip type operators. We are grateful to M. Haase for very useful com- 
ments on this theory, and refer to his book [27] for all the information necessary in 
this section. 

Let Ha = {z G C : \^z\ < a} denote the symmetric strip of height a > 0. If for 
an (unbounded) operator A on the Banach space X we have A G Strip(a), then we 
can apply the holomorphic functional calculus of A given by 

= dzf{z)D\{z,A),yF EMiUa], 

where -fa = dMa oriented positively (counterclockwise), and ^)l{z, A) is the resolvent 
of A for z G C. Define 

A{Ma) = {Fe Hol(Ha) : 37V G N s.t. F = 0{\^z\^)}, 

and 

b>a 

Now let I^Sj = — A + m*{t, x) be the known real lower semi-bounded operator 
acting on the vector bundle % over a Riemannian manifold Sj, and let C Et be 
a compact region. Denote 

V{K) = {/ G V{Tt) : supp/ C K}. 

Then we have the following result. 

Proposition Appendix C.l. For any F G -4[Ho] and f G 'D{K) it follows 

F{D^,)f eV{K). 

Proof. Let the nuclear topology be given by {Xi, 2, P), i.e., for any p G Nq we set 

(/, h)p = ^(g„,,(X,)/, QaAX^)h)L2 
q<p 

and consider the induced norms ||.||p. Define the Hilbert spaces 

then by the property of the countably normed spaces we have 

Hp C Hq, q <p, 

oo 

v{K) ^f]np. 

p=0 

Fix p, and define the operator Dp on Hp by setting Dpf = D^^f for all / G 
Dom(£)st)nHp, then Dom(Dp) D 'Hp+2 is a dense subspace of Hp. Then Dp is a real 
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symmetric operator, and hence by von Neumann's theorem possesses a self-adjoint 
extension Ap which needs not be lower semi-bounded. The self-adjoint operator Ap 

has a purely real spectrum, thus G Strip(O). Let y4,(Ha) 9 F{z) = 0(|3?z|^), then 
for a sufficiently large a < A € M, the function e{z) = {z — iA)"^^"'"^^ will regularize 
F on Ha. In particular, we will have [eF]{Ap) e B{Hp). Then F{Ap) = (Ap - 
a)^+2[eF](Ap) = [eF]{Ap){Ap-iX)^+'^, from where it follows that Dom(A^+2) c 
Dom{F{Ap)). From the definition of Ap it is clear that 'Hp+2{N+2) C Dom(A^+^), 
whence 'Hp^2(N+2) Goin{F{Ap)). Thus we have established, that whenever / e 
'Hp+2(N+2), then necessarily F{Ap)f e "Hp. Now if / € 'D{K), then for any p > 
we have / G 'Hp+2(N+2), and hence F{Ap)f E Hp. Meanwhile for any p > 0, the 
self-adjoint operator D^^ agrees with Ap on V{K). Therefore also their functional 
calcuH agree, F{Dj:Jf = F{Ap)f G Hp. Thus 

OO 

F{D,:,)fef]np=V{K), 
which completes the proof. □ 
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